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The self-field of a charged particle has a singular component that diverges at the particle. We 
use both coordinate and covariant approaches to compute an expansion of this singular field for 
particles in generic geodesic orbits about a Schwarzschild black hole for scalar, electromagnetic and 
gravitational cases. We check both approaches yield identical results and give, as an application, 
the calculation of previously unknown mode-sum regularisation parameters. 

In the so-called mode-sum regularization approach to self-force calculations, each mode of the 
retarded field is finite, while their sum diverges. The sum may be rendered finite and convergent 
by the subtraction of appropriate regularization parameters. Higher order parameters lead to faster 
convergence in the mode-sum. To demonstrate the significant benefit they yield, we use our newly 
derived parameters to calculate a highly accurate value of F r = 0.000013784482575667959(3) for 
the self-force on a scalar particle in a circular orbit around a Schwarzschild black hole. 

Finally, as a second example application of our high-order expansions, we compute high-order 
expressions for use in the effective source approach to self- force calculations. 



I. INTRODUCTION 



The notion of a self-force has a long history in physics. Early work led to the derivation of the Abraham-Lorentz- 
Dirac lj formula for the radiation reaction force on an accelerating point electric charge moving in a fiat spacetime. 
In the 1960s, DeWitt and Brehme [2] derived the curved spacetime equivalent, and a correction was later provided 
by Hobbs [3]. For a comprehensive review of the self-force problem, see Refs. [4j-[6]. 

With the advent of gravitational wave astronomy, the past two decades have seen a surge in interest in the self-force 
problem. This has been motivated by the study of so-called Extreme Mass Ratio Inspiral (EMRI) systems. These 
binary systems - consisting of a compact object of ~ 1 solar mass spiralling into a massive black hole of ~ 10 6 solar 
masses - are one of the most proimising candidates for study by planned space-based gravitational wave detectors 
[TM13] . The nature of interferometric gravitational wave detectors is that unlike traditional electromagnetic spectrum 
telescopes they are not directional; instead they see all directions at once. As a result, a crucial component of any 
observation is the use of data analysis techniques such as matched filtering to disentangle the desired signal from 
the noise. This, in turn, requires an accurate model of the signal one expects to see. This is even more true given 
the recent proposals for the eLISA/NGO detector P3| (which will have a reduced sensitivity compared to LISA, the 
previously proposed space-based gravitational wave detector) will mean that accurate waveform models are all the 
more crucial. 

The past two decades have seen new derivations of the equations of motion of a point charge moving in a curved 
spacetime; this was first done in the gravitational charge case by Mino, Sasaki and Tanaka [15] and Quinn and 
Wald 16 and in the minimally coupled scalar charge case by Quinn |17j . Since then, there have been significant 
achievements in putting these derivations on a firmer footing. The use of distributional sources in Einsteins equations 
is known to be problematic [TH]. Nevertheless, with sufficient care they have proven to be a useful practical tool 
in many cases [THl [20] ■ In the case of the self- force problem recent derivations have avoided the introduction of 
distributional sources altogether through the use of matched asymptotic expansions and careful limiting procedures 
|21H26j . At first perturbative order it is satisfying that these more rigorous derivations reproduce the same equations 
of motion as one would obtain using distributional sources. Nevertheless, it is likely that these more rigorous methods 
are required to advance to second perturbative order [2"7H3"T] and beyond [3"2"tI3"4"] . 
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Several practical self-force computation strategies have developed from these formal derivations, all of which are 
based on the now-justified assumption that the use of a distributional source is acceptable at first perturbative order: 

• The mode-sum approach [351 136j 

• The effective source approach 37, 38 

• The matched expansion approach [391 140] 

The key to all three approaches is the subtraction of an appropriate singular component (which only appears due to 
the use of a distributional source) from the retarded field to leave a finite regular field which is solely responsible for 
the self-force. This singular component must satisfy two criteria: 

1. It has the same singular structure as the full retarded field on the particle's worldline. 

2. It does not contribute to the self-force. 



There are many choices for a singular field that satisfies the above criteria, although not all choices are equal. Detweiler 
and Whiting [3T] identified a particularly appropriate choice. Through a Green function decomposition, they derived a 
singular field which not only satisfies the above two criteria, but also has the property that when it is subtracted from 
the full retarted field, it leaves a regularized field which is a solution to the homogeneous wave equation. Extensions of 
this idea to second perturbative order [29j [31] and fully non-perturbative contexts [34] have recently been developed. 

In this paper, we develop highly accurate approximations to the Detweiler- Whiting singular field of point scalar 
and electromagnetic charges as well as that of a point mass. This is achieved through high-order series expansions 
in a parameter e, which acts as a measure of distance from the particle's world-line. As examples, we derive explicit 
expressions for the case of geodesic motion in Schwarzschild spacetime and show how these may be applied to both the 
mode-sum and effective source approaches described above. Similar expansions may also be used for the quasi-local 
component in the matched expansion approach [42 44 whereby the Green function is matched onto a quasi-normal 
mode sum and the retarded field is then computed as the integral of this matched retarded Green function along the 
worldline of the particle. 

While the primary focus of this paper is on computing the singular field for Schwarzschild spacetime, many of the 
expressions we give are valid in more general spacetimes. In particular, where space allows, we do not make any 
assumptions about the spacetime being Ricci-flat. To make this distinction explicit, we use the Weyl tensor, C a bcd, 
in expressions which are valid only in vacuum and the Riemann tensor, R a bcd m expressions which are also valid for 
non-vacuum spacetimes. Note that this is done only for space reasons; our raw calculations include all non-vacuum 
terms in addition to those given in this paper and we have made the full expressions available in electronic form 45j. 

The layout of this paper is as follows. In S ec. [IT] we give exact formal expressions for the singular field in scalar, elec- 
tromagnetic and gravitational cases. In Sec. |III| we describe covariant and coordinate approaches to the computation 



of expansions of the key fundamental bitensors appearing in the formal expressions for the singular field. In Sec. IV 
we give explicit expressions for the singular field in the form of covariant series expansions. Our equivalent coordinate 
expressions are too large to be of use in print; instead we have made them available electronically [35] ■ In Sec. |VJ 
we use our coordinate expansions to derive high-order regularization parameters for use in the mode-sum method. In 
doing so, we give the next two previously unknown non-zero regularization parameters in scalar, electromagnetic and 
gravitational cases. In Sec. |VI| we apply our coordinate expansions to the effective source method and compute an 
effective source for our high order singular field. This improves on previous effective source calculations by adding an 



additional four orders, bringing the approximation to cight-from-leading order. In Sec. VII we summarize our results 
and discuss further prospects for their application. Finally, in Appendices [A] and [B] we give explicit expressions for 
many of the expansions in coordinate and covariant form, respectively, along with a more in-depth derivation of our 
expressions. 

Throughout this paper, we use units in which G = c = 1 and adopt the sign conventions of 46J. We denote 
symmetrization of indices using parenthesis (e.g. (ab)), anti-symmetrization using square brackets (e.g. [ab]), and 
exclude indices from (anti-)symmetrization by surrounding them by vertical bars (e.g. (a|6|c), [a|6|c]). We denote 
pairwise (anti-)symmetrization using an overbar, e.g. R^-^h — \{Rabcd + Rcdab)- Capital letters are used to denote 
the spinorial/tensorial indices appropriate to the field being considered. 

In many of our calculations, we have several spacetime points to be considered. Our convention is that (see Fig. [TJ: 

• the point x refers to the point where the field is evaluated 

• the point x' refers to an arbitrary spacetime point 

• the point x refers to an arbitrary point on the world-line 
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• the point X( a dv) refers to the advanced point of x on the world-line 

• the point Xr Iet \ refers to the retarded point of x on the world-line 

In computing expansions, we use e as an expansion parameter to denote the fundamental scale of separation, so that 
= 0(e). Where tensors are to be evaluated at one of these points, we decorate their indices 
appropriately using ' and" e.g. T a , T a and T a refer to tensors at x, x' and x, respectively. 



II. SINGULAR FIELD AND SELF-FORCE OF A POINT PARTICLE 



In an appropriate gauge, the retarded field, ip A (x), of an arbitrary point particle satisfies the inhomogeneous wave 
equation with a distributional source, 



2>W 



4ttQ / u a S 4 (x,z{t'))(It', 



where 



V A B = 5 A B {U - m 2 ) 



P A B 



(2.1) 



(2.2) 



is the wave operator, □ = g a6 V a Vb, g ah is the (contravariant) metric tensor, V a is the covariant derivative defined by 
a connection A A Ba- V a (^ A = d a ip A + A A Baf B , m is the mass of the field, Q is the charge of the particle, P A b(x) is 
a potential term and u A (x) is an appropriate tensor constructed from products of the four- velocity on the world- line 
parallel transported to x. The retarded solutions to this equation give rise to a field which exerts a self-force 



pa 



(2.3) 



on the particle, where p a a(x) is a tensor at x and depends on the type of charge. The distributional nature of the 
source leads to this self-force being divergent at the location of the particle and a regularisation scheme must be 
employed. We require a singular field, <f A S y that captures the singular behaviour of ip A and that, when subtracted 

from ip A , leaves a finite physical self- force. 

Detweiler and Whiting [41] showed how such a singular field can be constructed through a Green function decompo- 
sition. In four spacetime dimensions and within a normal neighborhood, the Green function for the retarded/advanced 
solutions to Eq. (2.1 ) may be given in Hadamard form, 



G( rct )/(adv) A B' (x, x') = 0_ /+ (x, x') {U A B > (s, x') 5 [a [x, x 1 )] - V A B < (x, x') 6 [-a (x, x')] } 



(2.4) 



where S [a (x, x')] is the covariant form of the Dirac delta function, 9 [a (x, x')] is the Heaviside step function, while 



U b' {x,x') and V b' (x,x') are symmetric bi-spinors/tensors which are regular for x' 



The bi-scalar a (x, x') 



is the Synge [4| world function; it is equal to one half of the squared geodesic distance between x and x'. The first 
term here, involving U A b' {x,x'), represents the direct part of the Green function while the second term, involving 
V A b' (x, x'), is known as the tail part of the Green function. 

Detweiler and Whiting proposed to define a singular Green function by taking the symmetric Green function, 
G( sym ) A B> = \{G(-[ C t) A B' + G( a dv) j4 B') and adding V A b> (x,x') (a homogeneous solution to Eq. (2.1 )). This leads to 
the singular Green function, 



G {S) A B > (x, x') = i { U A b' (x, x') 5 [a (x, x')] + V A B - {x, x') 6 [a (x, x')] } 



(2.5) 



Note that this has support on and outside the past and future light-cone (i.e. for points x and x' spatially separated) 
and is only uniquely defined provided x and x' are within a convex normal neighbourhood. Given this singular Green 
function, me may define the Detweiler- Whiting singular field, 



<P{S) 



T(rct) 



G (S) V {x,z{T'))u B 'dr', 



(2.6) 



T(adv) 



which also satisfies Eq. (2.1). Subtracting this singular field from the retarded field, we obtain the regularized field, 

V(R) = Vfret) ~ V(S)> ( 2 - 7 ) 
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which Detweiler and Whiting showed gives the correct finite physical self-force, 
Moreover, this regularized field is a solution of the homogeneous wave equation, 



v A Bvf K) = o. 



(2. 



(2.9) 



This holds independently of whether one is considering a scalar or electromagnetically charged point particle or a 
point mass. To make this more explicit, in the following subsections we give the form these expressions take in each 
of scalar, electromagnetic and gravitational cases. 



A. Scalar Case 

In the scalar case the singular field, $( s \ is a solution of the inhomogeneous scalar wave equation, 

(□ — £R — m 2 )$ (s) = q J V^Si(x, z{r))dT, (2.10) 

where q is the scalar charge and £ is the coupling to the background scalar curvature. An expression for <E>( S ) may be 
found by considering the scalar Green function [obtained by taking g A b 1 = 1 in Eq. (2.5)], 



G*( s ) = - {U{x, x')8[a(x, x')] + V(x, x')9[a{x, x')]} , 
which is a solution of the equation 

(D-ZR- m 2 )G (s) = -A%S(x, x'). 
Given this expression for the Green function, the scalar singular field is 

& S \x) =q f G {S) {x,z{ T ))dT 



U(x,x') 

<J c iU c ' 



+ 



V{x,z{r))di 



and one computes the scalar self-force from the regularized scalar field <I>( R ) = $( rot ) — $( s ) as 

pa _ g".b^) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



B. Electromagnetic Case 

In Lorenz gauge, the electromagnetic singular field satisfies the equation 



□ 4 S ) - i? a 6 4 S) = " 4 ^ e / 9aa'U a 'yf^6i(x,z(T))dT, 



(2.15) 



where e is the electric charge. An expression for A a may be found by considering the electromagnetic Green function 
[obtained by taking g A b 1 — g a b' in Eq. (2.5)], 



j {x, x 1 ) = - {U(x, x') aa ,S (a{x, x')) + V(x, x') aa ,9 {a(x, x'))} , 



where U a a ' = A 1 / 2 g a a < . This Green function is a solution of the equation 



(2.16) 



(2.17) 



Given this expression for the Green function, the electromagnetic singular field is 

J j 



u a U aa >(x,x') 



l-(adv) 



V aa ,{x,z(T))u a dr. 



T(ret) 



One computes the electromagnetic self- force from the electromagnetic regular field, — Aa Ct ^ — , as 

pa _ ab ca(R) 

9 A [ c ,b] ■ 



(2.18) 



(2.19) 



C. Gravitational Case 



In Lorenz gauge, the trace-reversed singular first order metric perturbation satisfies the equation 

OhfJ + 2C a c b d hf] = -167TA* / ga>( a u a ' g bW u b ' y/=g6t.{x, z{r))dT, (2.20) 



where fx is the mass of the particle and the trace-reversed singular field is related to the non-trace-reversed version 
by h-iT = ^ab ~ \h^g a b with = g ab h a S b ■ An expression for ft,^ ma Y be found by considering the gravitational 

(2.21) 



Green function [obtained by taking g A B' — g( a a'g b \' in Eq. (2.5)], 



G aba'b'( X i X ') = i {U(x,x')aba'b>5 [<t(x,x')] + V(x,x') aba ' b '0 [a(x, x')]} , 

with U ab a ' b i = A 1 / 2 g( a a ig b \i. This Green function is a solution of the equation 

OG ( a S b i b , + 2C a " b "G^ a , v = -4*g a , [agb) Mx,a/). 
Given this expression for the Green function, the trace-reversed singular first order metric perturbation is 



(2.22) 



h^=^ I G)Zb'(^z(T>))u a u b d T > 



•y(S) 



2fi 



u a u b U aba >v(x,x') 



T(adv) 



+ 2/1 / Vabafv (x, z{T))u a ' u b ' dr. 



(2.23) 



One computes the gravitational self-force from the regularized trace-reversed singular first order metric perturbation, 



ab ab 



pa = k abcd h m 



where 



6c;d' 



k abcd = \ ad u b u c _ ab u c u d _ I„a 4 o «I + + ^.fc 

2 2 4 4 



(2.24) 
(2.25) 



III. CO-VARIANT AND COORDINATE EXPANSION OF FUNDAMENTAL BITENSORS 

In the previous section, we gave expressions for the singular field in terms of the bitensors U A b 1 ( x , %') an d 
V A b' (x, x'). The first of these is given by 



U AB ' (x, x 1 ) = A 1 / 2 (a;, x 1 ) g AB ' (x, x') , 
where A(x,x') is the Van Vleck-Morette determinant [4], 



(3.1) 



A (x, x') = - [-g (x)}- 1/2 det (-a. afi , (x, x')) [-g (x')} 1/2 - det (- 9 Q ' Q (x, x') <r> (x, i')) , (3.2) 



6 



g AB is the bi-tensor of parallel transport appropriate to the tensorial nature of the field, eg. 

!1 (scalar) 
g ab ' (electromagnetic) (3.3) 

g a'(a g b)b' (g rav it a tional), 

and where the higher spin fields are taken in Lorentz gauge. Here, g a a (x,x') is the bi-vector of parallel transport 
defined by the transport equation 

<J a gab';a = = <J a 9ab';a'- (3.4) 

The bitensor V AB ' (x,x') may be expressed in terms of a formal expansion in increasing powers of a }47j : 

oo 

V AB ' (x, x') = J2 V » AB ' ^ x ') ° n s ') > ( 3 - 5 ) 
where the coefficients V AB (x,x') satisfy the recursion relations 

a ^ {A -l/2 v A B ' )ai + {n + 1) A -m V AB' + 1_ A -1/2 V B' ci yACl = ^ (3 _g a) 

for n G N, along with the 'initial condition' 

a ;c> (A -l/2 V AB'y al + A -1/2 V AB' + ^-1/2^^1/2^ = Q (3 6b) 

Looking at the above equations for U AB (x,x') and V n AB (x,x'), we see that a key component of the present 
work involves the computation of several fundamental bitensors, in particular, the world function a(x,x'), Van Vleck 
determinant A 1 / 2 (x, x'), four-velocity u a (x), and bivector of parallel transport g a b (x,x'). This may be achieved by 
expressing them as expansions about some arbitrary point x which is close to x and x' . We derive these here using 
both covariant and coordinate methods, each of which has its own advantages and disadvantages. The covariant 
expression is more elegant, allowing for compact formulas; however these formulas hide complex terms such as high 
order derivatives of the Weyl tensor that quickly become extremely time consuming to compute, even using computer 
tensor algebra packages such as GRTensorll 48 or xAct [13] ■ The coordinate approach is less elegant but more 
practical for explicit calculations and it avoids the need to use tensor algebra. Independently of the approach taken, 
these expansions may be used to compute expansions of U AB (x, x') and V n AB (x, x') (by substituting into the above 
equations), and hence of the singular field. In the case of covariant expansions, for explicit calculations one must 
further expand the covariant expressions in coordinates, yielding an expression which may be directly compared with 
those obtained from the coordinate approach. The resulting expressions are long but are explicit functions of the 
coordinates, enabling them to be transformed directly into, for example, C functions, indeed we give them in such 
form online |45j . 



A. Covariant Approach 



In this subsection, we briefly discuss our method for obtaining covariant expansions for the biscalars appearing in 
Eqs. (2.13), (2.18) and (2.23). We eventually seek expansions about a point x on the worldline (which we may treat 
as fixed in the majority of this paper). In doing so, we follow the strategy of Haas and Poisson |4l 150]: 

• For the generic biscalar A(x,z(t)), write it as A(r) = A(x,z(t)). 

• Compute the expansion about r = f. This takes the form 

A(t) = A(f) + A(f)(r -f)+ l -A{f){r - ff + • • ■ , 

where A(f) = A. a u a , A(f) = A. ah u a u h , 

• Compute the covariant expansions of the coefficients A(f), A(f), ■ ■ ■ about f. 



(3.7) 



7 




FIG. 1. We expand all bitensors (which are functions of x and x') about the arbitrary point x on the worldine. 

• Evaluate the expansion at the desired point, eg. A(x') = A(x, x'). 

• The resulting expansion depends on r through the powers of r — f. Replace these by their expansion in e (about 
x), the distance between x and the world- line. 

A key ingredient of this calculation is the expansion of At = r — f in e. The leading orders in this expansion were 
developed by Haas and Poisson [50] for the particular choices Ar + = v — f and Ar_ = u — r. They found 

Ar± = (f ± a) T ^ ^ Ruau* T [( f ± *)^— 1« - Ru*u*\*\ + 0(e 5 ), (3.8) 

where r = a a u a and s = (g ab + u a iP)<Ja<Ji. In Appendix [b] we extend their calculation to the higher orders required 
in the present work. In the same Appendix, we also apply the above method to compute covariant expansions of all 
quantities appearing in the expression for the singular field. 

In order to obtain explicit expressions, we substitute in the coordinate expansion for tr a (as discussed in Sec. [A]) 
along with the metric, Riemann tensor and 4- velocity (all evaluated at x). In doing so, we only have to keep terms 
that contribute up to the required order and truncate any higher order terms. 

B. Coordinate Approach 

In this subsection we follow a similar approach as described above, but in coordinates. We will start by considering 
two arbitrary points x and x' near x. We will seek expansion where the coefficients are evaluated at x, so we introduce 
the notation 

A/=/-i a , Sx a ' = x a ' - x a = x a ' - Ax a - x a , (3.9) 

where we use the convention that the index carries the information about the point: x a = x a . In the calculations 
below Ax a and 6x a are both assumed to be small, of order e. 

The first item we require for our calculations is a coordinate expansion of the biscalar a (x,x'), the Synge world 
function. We start with a standard coordinate series expansion about x, see for example |43j (note the difference in 
convention for Ax a in that paper), to get 

a(x, x') = y ab {x)Sx a ' Sx b ' + A abc {x)5x a ' 6x b ' 5x c ' + B abcd (x)Sx a ' 5x b ' Sx c ' 6x d ' +■■■ . (3.10) 

The coefficients arc readily determined in terms of derivatives of the metric at x by use of the defining identity 
2<T = o a <o a , see [43]. To be explicit, the first few are given by 

A a bc(x) = Z9(ab,c)(x) 

Babcd(x) = ^9(ab,cd){x) ~ j i 9 m {x) (g(ab.\p\{x)9cd),q{x) - ^9(ab,\p\{x)g\ q \c4){x) + ^g p (aA x )9\q\c,d)(x)) ■ 
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We now go one step further by expanding the coefficients about x to give a double expansion in Ax a and 5x a with 
coefficients at x. The first few terms are 

a(x,x') = \g- ai {x)8x a ' 8x b ' + \\g^{x)5x a ' 8x h ' Ax c + A- ais {x)8x a ' 8x b ' 8x c ' 

+ \z9ah,cd(x) 6xa ' Sxb ' ^ c Ax d + A- al - cd {x)8x a ' 8x b ' 8x c ' Ax d + B- al - cd {x) 8x a ' 8x b ' 8x c ' 8x d '] + 0(e 5 ), (3.11) 

where now we interpret 8x a as x a — Ax a — x a and we use square brackets to distinguish terms of different order in e. 
Rather than disturb the flow here and throughout this section we just give the first few terms of each expansion for a 
general metric to make the structure clear and give explicit expressions in Schwarzschild space-time to much higher 
order in Appendix |A"| 

Now that the coefficients are at the fixed point x, it is straightforward to take derivatives of cr at x and x' , for 
example, 



<V = 9ab$x b ' + [9ahJx b 'Ax c + 3A alE 8x b ' Sx c ' 

+ [y a - brcd Sx»'Ax c Ax d + 3A- ahe Jx»' 8x c ' Ax d + 4B- ahEd 8x b ' 8x c ' 8x d '] + 0(e 4 ), 
<J a = -<r a > + \g- b - c -Jx b ' 8x c ' + hg iEiS& 6x b ' 6x c ' 6x d ' + A- b5cl Jx b ' Sx c ' 8x d '} + <3(e 4 ), 



&a'b = -9ab ~ 9ab,c^ xC 



(3A aEd - 5 - 12B- ai5S )8x c '8x d ' - \g- alM Ax c Ax d + 0(e 3 ). 



Likewise we can calculate the Van Vleck-Morette determinant directly from its definition, 



As fx, a/) 



[-g(x)]-*\-v a , b (x 1 x')\[-g( X f )Y 



(3.12) 
(3.13) 
(3.14) 

(3.15) 



giving 



A 2 ( x , x') =1 + \g cd {2g azU - g a5g j - g^l) + \g cd g el (9ce,a9df,b + 2 gca,b{g S f,d ~ 2 9de,f) 
+ 29acA9bdj - 9bf,d) - 9ab,c(9ef,d ~ 2 9de,f)) &x a 8x b + 0(e 3 ). 



(3.16) 



To obtain an expressions at X( a d v ) an d ^(ret); we allow x a to be on the worldline and again give it as an expansion 
around the point x, as shown in Fig. [l] Writing x a in terms of proper time r gives 



x a (t) = x a + u a Ar + ±u a Ar 2 + ±u a Ar 3 



(3.17) 



where u a is the four velocity at the point x a and At = t — t. We are interested in determining the points on the 
worldline that are connected to a; by a null geodesic, that is we want to solve 

a(x a , x a '(T)) = =y a - b (u a Ar - Ax a ){u b At ~ Ax b ) 

+ [y a - b (u d Ar ~ Ax a )u b At 2 + \g- ab - c {u d Ar - Ax a )(u h Ar - Ax»)Ax c 

+ j9ab,M(u a ^r - Ax a )(u b AT - Ax b )(u 5 Ar - Ax c )] + 0(e 4 ). (3.18) 

By writing At = r±e + r 2 e 2 + r 3 e 3 + • • • and explicitly inserting an e in front of Ax a , we may equate coefficients of 
powers of e to obtain 

n 2 + 2g- al u d Ax b T 1 - g- ab Ax a Ax b = 0, (3.19) 
<fe E u s (u s Ti - Ax a )r 2 = - [iffas^n - Ax a )u b n 2 + i^f/n - Ax^uK, - Ax b )Ax c 

+ \g-ai Mi^n - Ax a )(u b n - Ax h )(u s n ~ Ax c )] . (3.20) 



Equation (3.19) is a quadratic with two real roots of opposite sign (for x spacelike separated from x) corresponding 
to the first approximation to our points X( a d v ) an d x (ret)i 



n± = g-aT,u a ^x b ± J{g- ab v?Ax b ) 2 + g- ab Ax a Ax b = f (1) ± p, 



(3.21) 



where rn\ is the leading order term in the coordinate expansion of the quantity r appearing in our covariant expansions. 
Equation (3.20) is typical of the higher order equations giving r n in terms of lower order terms. 
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IV. EXPANSIONS OF THE SINGULAR FIELD 

In this section we list the covariant form of the singular field to order e 4 , where e is the fundamental scale of 
separation, so, for example, f, s and a(x,x) a are all of leading order e. The cooordinate forms of these expansions 
are too long to be useful in print form so instead they are available to download [451. 



A. Scalar singular field 



To 0(e 4 



$(S) = Jl 

I s 



where 



the scalar singular field is 
1 



$(S)" 


= 15 


-2 -2 

r — s 




(3) 





6s 3 Gu ° u ° ' 24s 3 



Cu.(TU(tC< 



(r 3s )rC UJUI7;u -(r —5 )C U 



360s 5 



$(S) 



i 



(3) 4320s 5 



$(S)1 +0(e 5 )j, 

J (4) J 

(4.1) 



uaua^uaua 



(t s ) (SC^ct^ct^ctct -4- 4C li0 - - a C li(J(T ) -4- (r 6r s 3s )(4C U(Ttia C UCTU 



and 

$(S) 



(4) 



) \ +s \2C u a u [(r +s )C aSi ab 



+2f (f 2 + 3s*)C u - aa - b ] + 2C u \ b [2fC a - aa - b + (f 2 + s 2 )C ubaa ] + [f 4 + 6f 2 s" 2 + t]C u - aui C u a J> 
+ [r 2 + s 2 ]2C u&(J iC u a a b + C aacr i ) C a a c 



(4.2) 



30C* 



uaua 
2 

4 



f(3f 4 - lOfV + \bs A )c uaua , u - m{f z - -s z yc uaua , 

2l/~< or^2 -2i 

2 ^2\ 



•2s' 



-(- 3r[r' 3s ]Cuaua;uaa 3[r' S ]C uaua -.a a a 3[r 6r S 3s jC^o-^cj-^ua- 9?™[r lOr s 

^2 o-2\ 



3Cuaua;uuu r [ r ~ 10 J* S — 15s ] 
4 in^2-2 



15s \Cuaua\u C U aa [^^Cuaaa^a^ S ) T*(7* 3s ) (lOC^^^a^ loC^o-ua^ ^C U aua;a) 

-30C„ CTu5;u (f 4 - 6f 2 s 2 - 3s 4 )] - C uau a [S6r(r 2 ~ 3s 2 )C U aaa-a + (r 4 - 6f 2 s 2 - 3s 4 ) (13C uffua;a + 5C, 



u truer: a 



25C ' U(J (j a;u)] 



12s 4 



Co- cr [C^ag-^.g- T (Q oaob-.u 2C uao . fc;(T ) (r + s )(2C uaab . u C uau;);0 .) 
— f(f 2 + 3s 2 )C uatl 5 ;u ] + 2C tl a cr h [C . a(T 5 ;(T f + (f 2 + s 2 )(C uacr 5 ;0 . + C, u t ffa;[J — C CTacr 5 ;u ) 
+f(f 2 + 3s 2 )(C uati 5 ;cr — C u&tJ i. u — C ubaa . u ) — C^ui^f 4 + 6r 2 s 2 + s 4 )] + C tl a u b [(f 2 + s 2 )C . a(T b. (T 

— f(f 2 + 3s 2 )(C (jacr t; U — ^uaab;a) + + 6r 2 S 2 + S 4 ) (C U aub\a — 2C uao -t: tI ) 

-C uau ^ 4 + 10f 2 s 2 + 5s 4 )]l. (4.3) 



B. Electromagnetic singular field 



To C(e 4 ) 



the electromagnetic singular field is 



3vS Cauua H~ C U( j U( j{t S )w a 



a f u & , 1 

eja" — + 

I S DS- 

-\- ^t[t 3s *)Cuaua\u (j* ^ ) Cuauaa ^ 



+ 



24s 3 
1 

2880s 5 



4s (Cafi^fj^ ^Cauuu^fj^) 



A 



(S) 



1 



(3) 25920s 5 



.4 



(S) 



(4) 



O 



(£ 5 )}, (4.4) 



where 

[4 S) : 



(3) — ^20C u <jua;ua s {r + s 2 ) — 120C ucruaiaa rs + 120C n c M C a5aa >rs e — 2AQC uaua C U aua r s 2 {r 2 — 3s 2 ) 
- 360C UffaE C„ CT „ e s 4 (f 2 + s 2 ) - 120C tlffua; „ tl fs 4 (f 2 + 3s 2 ) + l20s 6 Cj/(C ucl - a - c f + C ad - aS ) 
+ 40C uEac[ Cu 5 u d s 6 (3f 2 + s 2 ) - ^OC^C^rs 4 + 2C uffu e s 4 (r 2 + s 2 )] - 120C uauE [C UCTff E s 4 (f 2 + s 2 ) 



2 (f 2 -3s 2 ) + 24C UCT 

+ 120C llCTilff C U(TUCT (f 2 - s 2 ) 2 + 16C u e JC aEaii s A (f 2 + s 2 ) + 2AC U aua;uuS 2 {f A - 6f 2 s 2 - 3s 4 

32C- n 5 ~ 2 ^ 4 C- 2 « 2 '5- 4 ^ i Qr< -n 5 (i-4^-4 , «-2-2 i -4^ , ,^-4^ c d 



„.„cC UCTU c ^(f 4 - 6f 2 s 2 - 3s 4 ) + 8C uSu(I Cu C u d s\f 4 + 6f 2 s 2 + s 4 ) + 16s 4 C u c /[2C (j5ff(i -f 
+ C uSa5 (f 2 + s 2 )] + 16s 4 C uSacI [Cu 5 J(f 2 + s 2 ) + 2C u 5 Jf(r 2 + 3s 2 )]K (4.5) 
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and 

)( s )l,., - 



[Aq ](4) — 216C U(TU a;a-(T(T rs H~ §^QC U(TU( j^(jC uaua rs (r 3s ) + 720C uaua - a C U(TU(T rs (r 3s ) 

rs (r + 3s 



ua\uua 

rs (r + 3s ) 

- 864C ucra c ;tl C tI(Ttie fs 4 (f 2 + 3s 2 ) — 216C U(Tua - uuu s 4 (f 4 + 6f 2 s 2 + s 4 ) — 540C tIcrttcr;u C t ( (TMa s 2 (f 4 
_ 6f 2 5 2 _ 3 - 4) + 720C UCTua;u C MffMCT s 2 (~f 4 + 6f 2 s 2 + 3s 4 ) - A32s 6 C a s a d , a (C u - cud f + C uda5 ) 

- 648s 6 C U C cr d :cr{C udas f + C adae ) + l6s S C U a Cd:e (C uS d g f — C uSsd f + C a£ds ~ C aSsd ) 

+ 432s 4 C uaa s , a [C uaa5 f + C U auc(r 2 + s 2 )] + 504s 4 C„ SCT % [C ua(T cf + 2C uau5 (f 2 + s 2 )] 

- 324s 4 C u ^ s . u [C u - at75 (r 2 + + C u - au - c f{f 2 + 3s 2 )] + l08s 4 C u ^. a [5C u - aaS (f 2 + s 2 ) + 6C UCTaE (f 2 + s 2 ) 

+ iC u - a u- c r{f 2 + 3s 2 )] - 96s 4 C uaa s . u [3C ua<j£ (f 2 + s 2 ) + 8C uiyu - c ?{f + 3s 2 )] + ms 4 C uau - a rc [iC uaiy - c {f 2 + s 2 ) 

+ 4C„ CT „ c -f(f 2 + 3s 2 )] + 24s 4 C ua „ E ;CT [9C U(TCT5 (f 2 + s 2 ) + 17C„™ E f(f 2 + 3s 2 )] - 216s 4 C uau s . u [3C uarT - c f{f 2 

+ 3s 2 ) + 6C u<T acf(r 2 + 3s 2 ) + 2C uauE s 4 (f 4 + 6f 2 s 2 + s 4 )] - 72s 4 C u - au 5 . u [8C uaa - c r (f 2 + 3s 2 ) 

+ 7C UffuE (f 4 + 6f 2 s 2 + s 4 )] - 216s 6 C u e /^[3C CTgad -f + C MEad -(3f 2 + s 2 )] - 144s 6 C ti E a ^[3C ttaVe f 

+ C uEud (3f 2 + s- 2 )] + 485 6 C ua< * d [3C u5r7d f + 3C uda - c f + 3C a5ad + C u - cud {3f 2 + s 2 )} 

+ 216s- 6 a c V; M [3C„ dCTE f + C u - cud (3f 2 + s 2 )] + 144s 6 C u V ; „[3C CTd - ae f + C ud - ae (3f 2 + s 2 )] 

+ 48s 6 C„ au E;J [3fC ffeff(i - + 3fC u5ud (f 2 + s 2 ) + C u - cad {3f 2 + s 2 ) + C uda5 (3f 2 + s 2 )} 

+ 2m 6 C u 5 a d ;u [3C u5ud f(f 2 + s 2 ) + C uda5 (3f 2 + s 2 )] + 144s 6 C^ ;u [3C„ EaJ f(f 2 + s 2 ) + C aEaa (3f 2 + s 2 )] 
+ {45C EjE ^C ed - s/ s 8 - 45C E ^„C ed - g/ 

>$QC U( j U( j : (j G a S (v ~ S ) H~ 36 Cu(7U(j ; uuuTS 

L (f 4 - 10f 2 s 2 - 15s 4 ) - 36C MffMCT;miff s 2 (f 4 - 6f 2 s 
+ 180C U(J11(7 ; 1i C 1i<71tff (3f 6 - 10f 3 s 2 + 15fs 4 ) - 216s 2 C ur7a 5 [rT [2fC urTu5 (f 2 - 3s 2 ) - 2C UCTffC -(-f 2 + s 2 )] 

- 72s- 4 C, E (T J ;CT [2C ug(T(i -r-+ C ffard - + C u5ud (f 2 + s 2 )} - lUs 4 C u s J, a [C a5rTd f+C u ^ d (f 2 + s 2 ) 

+ C uda5 (r 2 + s 2 ) + C uEud f(f 2 + 3s" 2 )] + 12s 4 C a \ d , u [C a - md f + 2C uSad (r 2 + s 2 ) + C u5ud f(f 2 + 3s 2 )} 
+ 60s 2 C uaa s . u [2C ua<T£ f{f 2 - 3s 2 ) + 5C uau£ (f 4 - 6f 2 s 2 - 3s 4 )] + 72s 2 C UffU E ;u [3C MCT „ g f(f 4 - 10f 2 s 2 

- 15s 4 ) + hC uaa - c (f - 6f 2 s 2 - 3s 4 )] - 72s 4 C u 5 J, a [C^ d (f 2 + s 2 ) + 2C u - cad f{r 2 + 3s 2 ) 

+ C ueud {f 4 + 6f 2 s 2 + s 4 )] + 144s 4 C M E CT J ;u [C CTECTd -(f 2 + s 2 ) + C uEad r(f 2 + 3s 2 ) + C udaE f(f 2 + 3s 2 ) 
+ C u5ud (f 4 + 6f 2 s 2 + s 4 )] - 60s 2 CU„ ff ;E [C UCTCTE f(f 2 - 3s 2 ) + C uauE (f 4 - 6f 2 s 2 - 3s 4 )] 

- 12s 2 C u ™V[16C„ CTCTC -f(f 2 - 3s 2 ) + 13C UCT „e(r 4 - 6f 2 s 2 - 3s" 4 )] + 12s 4 C u \ d , u [C a - cad f(f + 3s 2 ) 

+ 2C u5ad (f 4 + 6f 2 s 2 + s 4 ) + C u5ud f(f 4 + 10f 2 s 2 + 5s 4 )]}u a . (4.6) 



^ _ 

^2-2 



2\2 



3s 4 ) 



C. Gravitational singular field 

To (9(e ), the gravitational singular field is 



/i (S) 



a 9b \ — 



, '*s 1 

s 6s 3 



24s 3 



12s" 4 (C a 



(r 2 - s 2 )C uiyui yUaU- b - 6rs 2 C uau(a u- b) - 6s 4 C aubu + 
+ 4s 2 (u { aC- b)uua . u (f 2 + s 2 ) - fu { aC b)uua .Sj + u- a u b {r{f 2 - 3s 2 )C„ CTtlCT;M - (f - s 2 )C uaua -^ 

0(e 5 )}, 



aubu\<7 ^aubuiuJ 



l 



1440s 5 



ab 



+ 



l 



J (3) 6480s 5 



m 

ab 



J (4) 



(4.7) 



where 



(3) 



= -240C au ^ CT s 6 + 240C attl)ti;U(T fs 6 + 480C a 5 ua C buuS fs 6 - 120C aua s C buaE s 6 + 960C aub 5 C uau5 fs e 



+40C a „ 6u ' c E s 8 + 20C au cd C buSd s 8 + 240C' a c b d C uSud s 8 + 360C auua C buucr s 4 (r 2 + s 2 ) 
+240C aubu C uaur7 s 4 (f 2 + s 2 ) - 80C aubu , uu s 6 (3f 2 + s 2 ) - 40C a „ u E s 6 [6C 5 „ CTE f + C buu - C (3f 2 + s 2 )] 
+{120C^ ff;(TCT fs 4 + 2mC buua C uouo ft(f - 3s 2 ) - 120C 5 „ UCT;MCT s 4 (f 2 + s 2 ) - 360C 5 E MCT C UCTuE s 4 (f 2 + s 2 ) 
+l20C buua , uu fs 4 {f 2 + 3s 2 ) + l2QC- b % d [C u - cud f + C mEctJ ]s 6 + 120C 5 „/[C„ CT(7C -f + 2C uau5 (f 2 + s" 2 )]s 4 
+40C b s J[3C u5ad f + C usud (3f 2 + s 2 )]s 6 + mC^CU^f 2 + s 2 ) + C u „ uS f(f 2 + 3s 2 )]s 4 }u a 
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+{4C crgo . a C (T c ( j' i s 4 — 5C 5dS fC cde t s s + 12C M(TM( j ;(T( jS 2 (r 2 — s 2 ) + 16C llo - (J gC ucr(T c s 2 (r 2 — s 2 ) 
-12C„™ CT;tlff fs 2 (f 2 - 3s 2 ) + 64C u „ u B C ua „- c fs 2 (f 2 - 3s 2 ) + mC 2 uaua {f 2 - s 2 ) 2 + 8C/ u ^/(f 2 + s 2 ) 
+l2C uaua , uu s 2 {f i - 6f 2 s 2 - 3s 4 ) + 16C uau ,C uau 5 s 2 (f 4 - 6f 2 s 2 - 3s 4 ) + AC uEud C u 5 Js 4 (f 4 + 6f 2 s" 2 + s 4 ) 
+8C u e /[2C CTg ^f + C uda5 (f 2 + s 2 )]s 4 + 8C u5ad [C u s J(f 2 + s 2 ) + 2C u e M d >(F 2 + 3s 2 )]s 4 K^ (4.8) 



and 



m 

ab 



— 27Q{C aubu - aaa — C aubu;ucra r)s & — 1080(C ucru c - a C aube + C a c ua - a C buuS )rs 6 — 2700C atlf) c ;(T rC M(JMC s 6 

(4) 

— 90C aubu ,c Sl7 s s + OOCg^'^fs 8 — 360C U C U ^CW^s 8 + 720C U C U ;u fC a5bd s 8 + 180C u c a ; d C aube s 8 

+ l80C u C u d - 5 rC aubd S 8 ~ 120C au cd - a C buBd S 8 — 60C aua C ' d C buSd S 8 + 60C atl Cd ;urC bu5d S~ 8 — 180Ca C a d ; d C buuS S 8 

— l80C a c u d - d fC buuS s 8 — 720C a c b d - a C uSud s 8 + 360C a c ;, d -,urC uSud s 8 — 270C uau<T; crC aubu s 4 (f 2 + s 2 ) 
-1080C amff;CT C 6uuCT s 4 (f 2 + s 2 ) - 540C au5u;CT C UCTUff s 4 (f 2 + s 2 ) - 270C au6 „ ;tmu fs 6 (f 2 + s 2 ) 

+270(C uaua - u C aubu + 4:Cauua-;uCbuu<j + 2C aubu;M C ucrU(T )fS 4 (f 2 + 3s 2 ) + 5A0C aua - C - a (C buu5 f + ^buac)^ 

+90(C aubu:uua + 6C tIcrtt c ;tt C au 5 5 + 2Ca C U a:uCbuuc + 6C au ;, C ; u C Mffll5 ) (3f 2 + s 2 )s 6 + 60C autt c,<i (3C'5 Sua r 
+3Cbco-d — CbtigJ^)^ 8 + 180C a „5 c,<i (C uett jf + C ugcra )s 8 — 180C aucr c ;u [3Ct M(TE r + C buuS (3f 2 + s 2 )]s 6 
+180C Su „V[-3fC 6c -„ CT + 3fC 5 „ CT5 + C 5wte (3f 2 + s" 2 )]s 6 + 90C aubu > 5 [3C uaa5 f + 2C uau5 (3f 2 + s 2 )]s 6 
+180C- auu s lu [-3C buu5 f(f 2 + s 2 ) + 3C beua (3f 2 + s 2 ) - ^ U(Tg (3f 2 + s 2 )]s 6 

+ { — 108C5 ttM(j;(T(Tcr fs 4 + 3QC u c a d ' e C beds s 8 + 36C u c u d ' e fC bsdg s 8 + 12C b cde -aC u5ds s 8 + 24C b c rT d " e C u5ds s 8 
+l2C b cde . )U C u5ds fs 8 + 2AC b c u d ' e C uSdg fs 8 — 90(3C ucrua - a C buua + AC buua - a C uaU( j)rs 2 (f 2 — 3s 2 ) 
+108(C buua , uaa + 7C b 5 ua . a C uau£ )s 4 (f 2 + s 2 ) - 108(C buua . uua + 4C b 5 u ^, u C uau5 )fs 4 (f 2 + 3s 2 ) (4.9) 
+90(3C uaua - u C buurT + AC buua , u C uaua )s 2 {f 4 - 6f 2 s 2 - 3s 4 ) + lQ8C buua , uuu s 4 (f 4 + 6f 2 s 2 + s 4 ) 
—216C u c a d - cr (C bSud f + C bScrd )s 6 — 324C b c (7 d - (7 (C u5ud f + C uSad )s 6 + 8C bu cd ' e (C uSds f — C ugsd f 

+ C 'acde ~ C creed)* ~ 216C tt(Jcr C :<r[Ct M(T5 f + C buu5 (f + S )]s — 252Q mo . C ;( j [fC U(Jcre + 2C ucrttg (f +s)]s 

+MC urTU £ , a [6C b - cur7 (f 2 + s 2 ) - 5C bua5 (f 2 + s 2 ) - 3fC buuS (f 2 + 3s 2 )]s 4 + 162C uaa 5 ;u [C bua5 (f 2 + s 2 ) 
+C buu5 r(r 2 + 3s 2 )]s 4 + 12Q uu 5 ;ff [-9C MCTffe (f 2 + s 2 ) - lTfC^^f 2 + 3s 2 )]s 4 - A8C buua rc [3C ua<je {f 2 + s 2 ) 
+<±fC uau5 (f 2 + 3s 2 )]s 4 + 48C^ E ;u [3C Mffffe (f 2 + s 2 ) + 8C MCTUg f(f 2 + 3s 2 )]s 4 + 108C U<TU C , u [3fC bua - c {f 2 + 3s 2 ) 
-6C b5ua f(f 2 + 3s 2 ) + 2C buu5 (f 4 + 6f 2 s 2 + s 4 )]s 4 + mC buu 5 . u [8C uaa5 f{f 2 + 3s 2 ) + 7C uauE (f 4 + 6f 2 s 2 + s 4 )]s 4 
-72C u \ d , a [3C b5a<1 f + C b5ud -(3f 2 + s 2 )]s" 6 + 108Cj/. u [3C- b - cad r + C b5ua {3f 2 + s 2 )]s 6 
+108C u e u J ;u [3C 5EM(i -f(f 2 + s 2 ) + C bead -(3f 2 + s" 2 )]s 6 - 108^ e u J ;ff [3C Mgffd -f + C ugua (3f 2 + s 2 )]s 6 
+24C 5uCT c ^"[3C„ c - CTd -r- + 3C uSae f + 3C a - caS + C u5uS (3f 2 + s- 2 )]s" 6 + 72Cf a d , u [3fC uSaS + C u5uS (3f 2 + s 2 )]s 6 
+72C b \ d , u [iC ueuS f{r 2 + s- 2 ) + C u5ad (3f 2 + s 2 )]s e 

+2AC buu ^ d [3C a - caS f + 3C uSud f(f 2 + s 2 ) + C uead (3f 2 + s 2 ) + C udaS (3f 2 + s 2 )]s 6 }u- a 

+{ — 18(C (T CT -o- — C a a :«^)C 'o-co-js + 27C ^ ,jC adE js 8 — \8C~ ^ u^C sd ^jS 8 + 9C utTUtT - U(7rT fs (r — 3s ) 

— 15C M(TM(T ' fC uac! cS [f — 3s ) — 135C M(JM(T;(T C ucrtlcr (f — s ) J r9C uaua ^ uuu fs (f — lOr s — 15s ) 
+9C UffMff;utia s 2 (-f 4 + 6f 2 s 2 + 3s 4 ) + l5C u<Tua ' s C uau£ s 2 (-f 4 + 6f 2 s 2 + 3s 4 ) + C uaua , acja (-9f 2 s 2 + 9s 4 ) 
+45C MCTMCT;u C U(TU(T (3f 5 - 10f 3 s 2 + 15fs 4 ) - 54C uaa s . a [2C uaue r(f 2 - 3s 2 ) + C uaa - C {f 2 - s 2 )]s 2 
+36C MgCTd -s- 4 ha E a^f + C a \ d , u {f + s 2 )] - 18C Mgud -s 4 [a e /^(f 2 + s 2 ) - C CT E ff J ; „f(f 2 + 3s 2 )] 
-3&C u \ 3 , a [C a - c<yd r + C u - cad {f + s 2 ) + C u - da - c (f + s 2 ) + C u5ud f(f 2 + 3s" 2 )]s 4 + \hC ua S , u \2C uaa - c f (f 2 - 3s 2 ) 
+5C Ufft(g (f 4 - 6f 2 s 2 - 3s 4 )]s 2 + \8C UGU s , u [3fC uau - c {f 4 - 10f 2 s 2 - 15s 4 ) + 5C U(TCTg (f 4 - 6f 2 s 2 - 3s 4 )]s 2 
-\8C u \ d , a \C a - cad {f 2 + s 2 ) + 2C u - cad f(f 2 + 3s 2 ) + 18C ug „ a (f 4 + 6f 2 s 2 + s 4 )]s" 4 
+36C M E CT J ;u [C ffgad -(f 2 + s- 2 ) + C ugffa f(f 2 + 3s" 2 ) + C uSr75 f(f 2 + 3s" 2 ) + C u - cud (f 4 + 6f 2 s" 2 + s" 4 )]s" 4 
+C UCTU e :ff [-48C UCTffg f(f 2 - 3s 2 ) + 39CU„ g (-f 4 + 6f 2 s 2 + 3s 4 )]s 2 + 18C u 5 J ;u [C a5ad f(f 2 + 3s 2 ) 
+2C u - cad {f 4 + 6f 2 s" 2 + s- 4 ) + C u5ud f(f 4 + 10f 2 s" 2 + 5s" 4 )]s- 4 } WaU ,. (4.10) 
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V. MODE-SUM REGULARIZATION PARAMETERS 



The singular field expansions derived in the previous sections have several applications in explicit self-force calcu- 
lations. One of the most successful computational approaches to date is the mode-sum scheme of Barack and Ori 
[551 [35]; the majority of existing calculations are based on it in one form or another [5Url73| . The basic idea is to 
decompose the singular retarded field into spherical harmonic modes which are continuous and finite in general for the 
scalar case and in Lorenz gauge for the electromagnetic and gravitational cases. A key component of the calculation 
involves the subtraction of so-called regularization parameters - analytically derived expressions which render the 
formally divergent sum over spherical harmonic modes finite. In this section, we derive these parameters from our 
singular field expressions and show how they may be used to compute the self-force with extremely high accuracy. 



A. Rotated Coordinates 



In order to obtain expressions which are readily written as mode-sums, previous calculations [35j [50l [53] found it 
useful to work in a rotated coordinate frame. We found it most efficient to carry out this rotation prior to doing any 
calculations. To this end, we introduce Ricmann normal coordinates on the 2-sphere at x in the form 



W\ = 2 sin ^ — J cos/3, 
where a and /3 are rotated angular coordinates given by 



sin cos (, 
sin 9 sin <; 
cos ( 



w 2 = 2 sin ( — ) sin ,3, 



cos a, 

sin a cos /3, 

sin a sin (3. 



In these coordinates, the Schwarzschild metric is given by the line element 



ds 2 = - 



r — 2m 
r 



dt 2 



r 

r — 2m 



dr 2 



+ 2dw\dw2 



WiW 2 



4(4 - w\ - w 2 ) 



16 - w\ (8 - 

4 (4 - w\ 

16 - w\ (8 - 



4 (4 - wf - wl) 



dw\ 



dlVn 



(5.1) 



(5.2) 
(5.3) 
(5.4) 



(5.5) 



The algebraic form of the metric makes it very suitable for using with computer algebra programmes such as Mathe- 
matica. The apparent complexity of having a non-diagonal metric on S 2 is in fact minimal since the determinant of 
that metric is simply 1. 



B. Mode decomposition 

The method of regularization of the self force through Z-mode decomposition is by now standard, see, for example, 
[35], [53] and [50]. Having calculated the singular field, it is straightforward to calculate the singular self-forceQ F a , 
for scalar, electromagnetic and gravitational cases using Eqs. (2.14), (2.19) and (2.24) respectively. We study the 
multipole decomposition of F a by writing 



F a (r,t,a,(3) = J2F l a m (r,t)Y lm (a, (3) . 



(5.6) 



I in 



where Y 



lm , 



are scalar spherical harmonics, and accordingly 



FT (r, t) = / F a (r, t, a, 0) Y lm * (a, (3) dtt 



(5.7) 



1 In this section, for notational convenience we drop the implied (S) superscript denoting "singular" as we are always referring to the 
singular component. 
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To calculate the /-mode contribution at x = (to,ro,ao, Po), we have 

F l a (r , t Q ) = lim V F<™ (r + Ar, t„) F im («o, M • 



(5.8) 



In previous calculations, Eq. (5.6) has naturally arisen in Schwarzschild coordinates with 9q — ^ , and it was 
necessary to perform a rotation to move the coordinate location of the particle from the equatorial plane to a pole in 
the new coordinate system. However, by choosing to work in an S 2 Riemann normal coordinate system from the start, 
our particle is already located on the pole. This saves us from further transformation and expansions at this stage. 
With the particle on the pole, Y lm (a — 0,/3 ) = for all This also allows us, without loss of generality, to 

take fio — 0. Taking a- , /3q and m all to be equal to zero in Eq. (5.8) gives us 



(r + Ar,t ) 



21 + 1 



47T A 

For each spin field, the singular self-force, F a (r, t, a, f3), has the form 

„(3n-2) 

F a M, «,/?) = £ ^krr-^" 3 . 

n=l H 



lim F a (r + Ar,t , a, /3)P t (cos a) dQ. (5.9) 



(5.10) 



where = b aia2 ,,, ak (x)Ax ai Ax a2 . . . Ax ak . On identifying t\ — f(i) ± p, this form can be easily seen to follow from 
the coordinate representation of the above expressions for the singular field. In using Eq. (5.10) to determine the 
regularization parameters, we only need to take the sum to the appropriate order: n = 1 for A ai n = 2 for B a , etc. 
Explicitly, in our coordinates p = \J (g & iu a Ax b ) 2 + g a i,Ax a Ax b takes the form 



p(r,t,a,PY 



(E 2 r 3 Q -L 2 (r -2M)) 
r (r - 2M) 2 



Ar 2 



r o) A 



2Er f 
r a - 2M 



Ar + 2ELAwi I At 



2Lr r 
r a - 2M 



ArAw! 



2M 



E 2 + — 1 ) At 2 + r 2 Aw 2 2 , 



(5.11) 



where the a, /3 dependence is contained exclusively in Aw\ and Aw2- In particular, taking t — to (At — 0) allows us 
to write 



p(r, t Q ,a,PY 



We further define 



E 2 4 



(L 2 + r 2 )(r -2M} 



; Ar 2 + (L 2 + r 2 ) I A Wl 



Lr r 



(r -2M) (L 2 +r 2 ) 



Ar) + r 2 Aw 2 . (5.12) 



po (a, /3) 2 = p (r , t , a, P) 2 = (L 2 + r 2 ) Aw 2 + r 2 Ai 



and this allows us to rewrite our Aw's in the alternate forms 

2 



Awf = 2(1- cos a) cos (3 



Awl = 2(l-cosa)sin 2 /3 



Po 



(r 2 + L 2 ) X 
i 



Pa 



(r 2 + L 2 ) X 



2 

) 

2 

L 2 X 



cos 2 /3=^-(fc-(l- X )) 



sin 2 P= ^-(1 - X ), 



(5.13) 

(5.14) 
(5.15) 



where 



x(/3) = l-fcsin 2 /3, 



L 2 ' 



(5.16) 



Suppose, for the moment, that we may take the limit in Eq. (5.9) through the integral sign, then using our alternate 
forms we have 



am — „ , , e ^rr-, e — Po e c a(n )(r ,xj- 



Ar-yO pSn+l 



Pa 



2)1+1 



(5.17) 



14 



In [73], it was shown that the integral and limit in Eq.(5.9) are indeed interchangeable for all orders except the leading 
order, n = 1 term, where the limiting po~ 3 would not be integrable. Thus we find the singular self-force now has the 
form 



F l a (r ,t ,a,l3) 



21 + 1 

47T 



lim f Ei ^-P l (cos a) dn + V e"- 3 f p n ~ 3 c a(n) (r ,x) P 



(cos a) dfl 



= P l a [-i] (ro, to) £~ 2 + F l a[0] (r 0> to) e" 1 + F l a[2] (r 0> t ) e 1 + ^[4] (r , *o) e 3 + ^ [6] (r , *o) ^ + ■ ■ ■ , 

(5.18) 

where the f3 dependence in the c„'s are hidden in x, while the a and j3 dependence of F a (r, to, a, j$) is hidden in both 
the p's and c„'s. Note here that we use the convention that a subscript in square brackets denotes the term which 
will contribute at that order in 1/1. Furthermore the integrand in the summation is odd or even under Awi — > —Auii 
according to whether n (and so 3n — 2) is odd or even. This means only the even terms are non-vanishing, while 
^i[i](ro,to) = ^ [ 3 ] (ro ) to) = i ;, i [ B](ro,*o) = Oetc. 

Some care is required in dealing with taking the limit in the first term. As this has been addressed previously 
[5T71 1551 [741 [75] and our main interest is in the higher orders we omit the details here and only discuss the calculation 
of the higher order terms. 

In the higher order terns in Eq. ( |5.18 ), we may immediately work with po 2 = 2\(L 2 + rg)(l — cos a) so, 



Po 



(r ,td, a, ft" = [2 X (L 2 + r 2 ) (1 - cosa)]" /2 = [2 X (L 2 + r 2 )] n/2 ^^ /2 (0)fl (cos a) 



(5.19) 



where A t 2 (0) = \/2 from the generating function of the Legendre polynomials and, as derived in Appendix D of [53] , 
for(n + l)/2eN 



V n/2 = (_i)(»+D/ 2 2 1 +"/ 2 (n!!) 2 , 
/a w _ Pn/2 (21 + 1) 



A" /z (o) 



1 (2Z-ra)(2Z-n + 2)...(2Z + n)(2Z + n + 2) 

In this case the angular integrals involve 

V«/2 (/?) \ ^ 



1 



dp 



X. 



n 1 

-r~2^ k 



(5.20) 



(5.21) 



2W xW) n/2 

where (n + l)/2 G N U {0}. The resulting equations can then be tidied up using the following special cases of 
hypergeometric functions 

2 

— ; 

7T 

2 



^_i(fc)= 2 Fi 



2'2' i,fc 



2'2 ;1;fc 



: /C, 



(5.22) 
(5.23) 



where 



K 



7T/2 



(1 - fcsin 2 ^)- 1/2 (i/3, £ 



r/2 



(1 - fcsin 2 /3) 1/2 d/3 



(5.24) 



are complete elliptic integrals of the first and second kinds, respectively. All other powers of x can he integrated to 
hypergeometric functions that can then be manipulated to be one of the above by the use of the recurrence relation 
in Eq. (15.2.10) of 7(i . That is 



p - 1 r^ 1 W + 1 -^t (p - |) ^ W- 



P (k-iy 



p(*;-i) 



(5.25) 



In the next sections, we give the results of applying this calculation to each of scalar, electromagnetic and gravita- 
tional cases in turn. In doing so, we omit the explicit dependence on I which in each case is 



F. 



l [ _ 1] = (2l + l)F a[ . 1] , 



F' 



a[Q] 



F 



a[0]i 



F, 



,,[2] 



(21 - l)(2Z + 3)' 



F; 



F 



a[4] 



a[4] 



(2Z-3)(2Z- l)(2Z + 3)(2Z + 5)' 



F a[6] ~ 



F 



a[6] 



(21 - 5)(2Z - 3)(2Z - 1)(2Z + 3)(2Z + 5)(2Z + 7) ' 



(5.26) 
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It is also also worth pointing out that there exists in the literature several different notations for the regularization 
parameters. We have adopted a notation which is readily extensible to other orders and which makes the dependence 
on I explicit. To avoid confusion, in Table [T] we give the relation between our notation and other common notations. 



RP 


BO 


DMW 




A a 


A a 


F a [ ] 


B a 


B a 


Fall] 


C a 


C a 


F a {2] 




D a 


F a{4] 




El 


F a[6 ] 




El 



TABLE I. Relation between notational choices for the regularization parameters (RPs). The most common choices are those 
of either Barack and Ori [35 or Detweiler, Messaritaki and Whiting [53 . 



C. Scalar case 

In the scalar case, the regularization parameters are given by 

r sgn(Ar) _ Er sgn(Ar) 



(5.27) 



ft 



t[0] 



Er r 



;(2£-JC), 



F, 



n(L 2 + r 2 ) 3 / 2 

= 7rro(ro _ 2M i )(L2+ro2) 3 /2 {[ 2 ^ 3 - <*> - 2M ^ + r ^ £ - ^ + ^ - 2M ^ 2 + 

F m = 0, F m = - z -^^_ (£-K), 



(5.28) 



(5.29) 



where 



where 



F m = 



Er 

2nr 4 (L 2 + r 2 )V 2 



( F t[2]£ + F t[2]K)> 



(5.30) 



F t%] 



8E 2 {L 2 - rftrl - (L + r 2 ){36L b M + l04L 4 Mr 2 + 



98L 2 Mr 4 



I?rl 



F^ 2] = -E 2 r 7 {3L 2 - brl) + 2r 2 {L 2 + rg)(9i 4 M + 18L 2 Mr 2 + 13Mr% - 2r\ 



46Mr 6 -7r 7 ), 



0J> 



r[2] 



27rr6(r -2M)(L 2 + r2)7/2 



( F r[2]£ + F rl2]tyi 



r{2]> 



(5.31) 



F* [2] = —&E vq(L - 4) + AE z rl{L z + r^)(9L°M + 2QL 4 Mrf ) + 23L^Mr 4 

- (r - 2M){L 2 + rl) 2 (28L 6 M + 82L 4 Mr 2 + 82L 2 Mr 4 - L 2 r^ + 32Mr£ 



r2 5 



UMr 6 - 3r 7 ) 



F r \2] 



E 4 r^{2,L 2 - 5r 2 ) - E 2 rl(L 2 + r 2 )(18L 4 M + ML 2 Mr 2 + L 2 rf 



32Mr^ 



- 3r 7 ), 
7r 5 ) 



+ (r - 2M)rl{L z + rg) 2 (14L 4 M + 28L z Mr^ + WMr^ - rg), 



F 6[2] — 0- 



(5.32) 
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Fm " 2nLr*(J+r*)S/i {F *M £ + F *^ K) ' (5 - 33) 

where 

F%[2] = E 2 r 7 (7L 2 - r 2 Q ) + (L 2 + r 2 )(28L 6 M + 58L 4 Mr 2 + ML 2 Mr\ - L 2 r% + r 7 ), 
Fm = -E 2 rl(3L 2 r 2 ) r 2 (L 2 + r 2 )(UL*M + 16L 2 Mr 2 + r 5 ), 



where 



F5 4 , = -30^ 4 ri 6 (23L 4 - 82Z/V 2 + 23r 4 ) 



2E 2 r%{L 2 + r 2 )(44800L 12 M + 219136L 10 Mr 2 + 428252L 8 Mr% + 418776L 6 Mr® 
+ 206374L 4 Mr 8 + 45L 4 r 9 + 45188L 2 Mr 1 , - l230L 2 r^ - 166Mr 42 + 645r 43 ) 



- 2(L 2 + rg) 2 (20480L 14 M 2 - 9728(XL 12 M 2 r 2 , + 8512QL 12 Mr% - 700832L 10 M 2 r 4 
+ 388480L 10 Mr^ - 1426472L 8 M 2 r^ + 704552L 8 MrJ - 1358276L 6 M 2 r 8 1 
+ 635226L 6 M^ - 635180L 4 M 2 r^° + 286498L 4 Mr 1 , 1 - 15L 4 r^ 2 - 124540£ 2 M 2 r 1 l 2 
+ 54086L 2 Mr 1 , 3 - 90i 2 r^ 4 - 2796M 2 r 1 l 4 + 182Mr^ 5 + 285r^ 6 ), 

F*fa = 15£ 4 r 46 (15L 4 - 82L 2 r 2 + 31r 4 ) 

- AE 2 rl{L 2 + r 2 )(11200i 10 M + 44984L 8 Mr 2 + 68227L 6 Mr 4 

+ 46849L 4 Mr® + 13493L 2 Mr 8 - 270L 2 r 9 l + 127 Mr™ + 210r^) 
+ rl{L 2 + rg) 2 (20480L 12 M 2 - 115200L 10 M 2 r^ + 85120L 10 Mr% - 599072L 8 M 2 r\ 

+ 31400(XL 8 M^ - 908104L 6 M 2 r^ + 433792L 6 Mr£ - 589164L 4 M 2 rg + 268648L 4 Mt^ 
- 151484 J L 2 M 2 r 1 l ° + 66380L 2 Mr^ - 15L 2 ^ 2 - 5592M 2 r 1 , 2 + 1204M7- 1 , 3 + 345r 1 l 4 ), 



F rW ~ 40^ 3 (r -2M)(L2+r 2 )iV2^[4] f + F r[^)> (5- 

where 

F^ [A] = 30£; 6 r 49 (23L 4 - 82i 2 r 2 + 23r 4 ) 

- E A rl{L 2 + r 2 )(89600L 12 M + 438272L 10 Mr 2 + 856504L 8 Mr 4 + 837552L 6 Mr® 
+ 411938 J L 4 Mr 8 ) + 495L 4 ^ + 92836L 2 Mr 1 , - 3690i 2 r^ - 902Mr^ 2 + 1575^ 3 ) 

+ 8E 2 4(L 2 + rg) 2 (5120L 14 M 2 - 35200L 12 M 2 r 2 , + 2672(XL 12 Mr 3 1 - 227368L 10 M 2 r 4 + 123200L 10 Mr% 

- 456300L 8 M 2 r^ + 225979L 8 Mr 7 - 434510L 6 M 2 r 8 + 206277L 6 Mr 9 l - 203983 J L 4 M 2 r 1 ,° 

+ 94211L 4 Mr 41 - 40376L 2 M 2 r 42 + 18367L 2 Mr 1 , 3 - 135L 2 r 44 - 571M 2 r 1 1 4 - 146Mr^ 5 + 135r^ 6 ) 

- (r - 2M)(L 2 + r^) 3 (40960L 14 M 2 - 86016L 12 M 2 r 2 + 116480L 12 Mr^ - 860224L 10 M 2 r 4 
+ 510080L 10 Mr^ - 1780112L 8 M 2 r^ + 882400L 8 Mr\ - 1657392L 6 M 2 r 8 1 + 752340L 6 Mr% 

- 743164 J L 4 M 2 r 1 ,° + 316100L 4 Mr^ + 30i 4 r^ 2 - 136236L 2 M 2 r^ 2 + 53200L 2 Mr 1 , 3 + 75£ 2 r 1 1 4 

- 3120M 2 /- 1 , 4 + 160Mr 45 + 165r 46 ), 
= -15S 6 r 49 (15L 4 - 82L 2 r 2 + 31r 4 ) 

+ E 4 r™(L 2 + r 2 )(44800L 10 M + 179936i 8 Mr 2 + 272908L 6 Mr 4 

+187126L 4 Mr« + 135L 4 r£ + 55232L 2 Mr 8 - 1710L 2 r 9 + 118Mr 40 + 1035r 41 ) 

- E 2 r\{L 2 + r 2 ) 2 (20480L 12 M 2 - 158720L 10 M 2 r 2 + 106880L 10 Mr^ - 769632L 8 M 2 r 4 

+ 399280L 8 Mr% - 1159632L 6 M 2 r^ + 559556L 6 MrJ - 755876L 4 M 2 rg + 352004L 4 Mr^ 
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- 196524L 2 M 2 rl° + 89680L 2 MrJ 1 - 405L 2 r 1 1 2 - 5528M 2 rJ 2 + 512Mr^ 3 + 675rJ 4 ) 

+ (r - 2M)r 2 (L 2 + r^) 3 (20480L 12 M 2 - 60928£ 10 M 2 rg + 58240L 10 Mr 3 ) - 375840L 8 M 2 r 4 
+ 204080L 8 Mr§ - 564472i 6 M 2 r^ + 265360L 6 Mr 7 l - 350956L 4 M 2 r 8 + 152380Z 4 Mr 9 , 

- 84276L 2 M 2 r 1 1 ° + 33740L 2 Mr 1 , 1 + 15L 2 rJ 2 - 3120M 2 rJ 2 + 640MrJ 3 + 75rJ 4 ), 



F m = 0, (5.36) 



Ff [A] = -lS^rJ^SL 4 - 82L 2 r 2 + 3r 4 ) 

- 10£ 2 r^(£ 2 + r 2 )(4352L 12 M + 16512L 10 Mr 2 + 22948L 8 Mr 4 

+ 13346£ 6 Mr[> + 2136£ 4 Mr 8 - 9L 4 r 9 - 710L 2 MrJ" + 126L 2 r^ - 9rJ 3 ) 
+ (L 2 + r 2 ) 2 (40960L 14 M 2 - 96256L 12 M 2 r 2 + 116480L 12 M^ - 704064L 10 M 2 r 4 

+ 429440L 10 Mr^ - 1134992L 8 M 2 r^ + 595040L 8 MrJ - 755632L 6 M 2 r 8 , 

+ 372500L 6 Mrl - 194724L 4 M 2 r 1 ,° + 94940L 4 Mr 1 , 1 + 30£ 4 rJ 2 - 6276L 2 M 2 rJ 2 

+ 4040L 2 Mr ( 1 l 3 + 105L 2 rJ 4 + 480M 2 r 1 , 4 - 45rJ 6 ), 
^ 4] =15^ 6 (L 2 -3r 2 )(15i 2 -r 2 ) 

+ WE 2 r 7 (L 2 + r 2 )(2176L 10 M + 6352L 8 Mr 2 + 6018L 6 Mr 4 

+ 1666i 4 Mr^ - 320L 2 Mr 8 + 63L 2 r 9 - 9r^) 

- r 2 (L 2 + r 2 ) 2 (20480L 12 Af 2 - 66048L 10 M 2 r 2 + 58240L 10 Mrg - 293280L 8 M 2 r 4 
+ 163760 J L 8 Afr 5 ) - 314392L 6 M 2 rg + 156960L 6 MrJ - 114876L 4 M 2 r 8 l 

+ 55420L 4 Mr 9 , - 6516L 2 M 2 rJ° + 3740L 2 AlrJ 1 + 15£ 2 rJ 2 + 480M 2 r 1 , 2 - 45rJ 4 ), 



F ^ = , m J^l ^ 2 ( F ^ £ + F ^ K ) ' (5 - 38) 
56071-^ (L 2 + r^) v ' 



28000£ 6 r 23 (r - L)(L + r ) (llL 4 - 74L 2 r 2 + llr 4 ) 

- 25£ 4 r 8 (L 2 + r 2 ) ( - 16056320L 18 M - 10715136(XL 16 Mr 2 - 302586880L 14 Mr 4 - 464979968L 12 Mrg 

- 412568652i 10 Mrg - 201055024L 8 MrJ - 3926841(XL 6 MrJ 2 - 1575L 6 rJ 3 + 5226426L 4 MrJ 4 
+ 99435L 4 rJ 5 + 3185118L 2 AlrJ 6 - 186165L 2 rJ 7 + 19662MrJ 8 + 35385rJ 9 ) 

+ 2E 2 4 (L 2 + rl) 2 ( - 1007616000L 20 M 2 - 2885324800L 18 M 2 r 2 - 1548288000L 18 Mr 3 , 

+ 5271990272£ 16 M 2 r 4 - 9940582400L 16 Mr^ + 35832487264L 14 M 2 r^ - 27145052800L 14 MrJ 
+ 69571689904L 12 M 2 rg - 40793731200L 12 Mrg + 69887626312 J L 10 M 2 rJ° - 36329433800L 10 MrJ 1 
+ 39015325900L 8 M 2 rJ 2 - 19063343950L 8 MrJ 3 + 11166709052 J L 6 M 2 rJ 4 - 5373108900L 6 MrJ 5 
+ 7875L 6 rJ 6 + 1046817944L 4 M 2 rJ 6 - 575985300L 4 MrJ 7 + 94500L 4 rJ 8 - 118281276L 2 M 2 rJ 8 
+ 29440500L 2 MrJ 9 - 1178625L 2 r 20 - 8271468M 2 r 20 + 4798 50Mr 21 + 414750r 22 ) 

+ (L 2 +rlf ( - 5775360000L 20 M 3 + 2580480000L 20 M 2 r - 18980904960L 18 M 3 r 2 + 750796800L 18 M 2 r 3 1 
+ 342988800(XL 18 Mr 4 + 10876463104L 16 M 3 r 4 - 5306391552(XL 16 M 2 r 5 1 + 21396480000L 16 Alr<> 
+ 143196789568L 14 M 3 r^ - 191859546624L 14 M 2 r^ + 56793312800L 14 Mr 8 1 + 292841560608L 12 M 3 r 8 1 

- 317782413664L 12 M 2 rg + 83096000800L 12 MrJ + 297880915104 J L 10 M 3 rJ° - 295661821784L 10 M 2 rJ 1 
+ 72364880400 J L 10 MrJ 2 + 168534399040L 8 M 3 rJ 2 - 159472848000L 8 M 2 rJ 3 + 37560515600L 8 MrJ 4 
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+ 50707761864L 6 M 3 rJ 4 - 4682664082(XL 6 M 2 r 1 1 5 + 10839698800L 6 Mr^ 6 + 7000 J L 6 r 1 , 7 

+ 6272875728L 4 M 3 r 1 1 6 - 5878415984L 4 M 2 rJ 7 + 1398754050L 4 Mr^ 8 + 41125LV 1 , 9 - 86931352L 2 M 3 r 1 , 8 

+ 212436L 2 M 2 r 1 , 9 + 21357300L 2 Mr 20 + 124 2 50L 2 r 21 - 29620256M 3 r 20 + 16081176M 2 r 21 

- 597150Mr 22 - 24 5 8 75r 23 ), 

Fjfa = -875£ 6 r 23 (-105L 6 + 1189L 4 r 2 - 1531i 2 r 4 + 247rg) 

+ 50E 4 r^° (L 2 + r 2 ) ( - 4014080L 16 M - 23275520L 14 Mr 2 - 55468800L 12 Mr 4 - 68718512L 10 Mr^ 

- 45183275L 8 Mr 8 , - 13052460L 6 Mrj" + 362358L 4 MrJ 2 + 18900LV 1 , 3 + 919476L 2 Mr 1 , 4 - 49560L 2 r 1 , 5 
+15501Mr 1 , 6 + 12180r^ 7 ) 

- E 2 r% (L 2 + r 2 ) 2 ( - 1007616000L 18 M 2 - 200366080(XL 16 M 2 r 2 - 1548288000L 16 Mr^ 

+ 6977961472L 14 M 2 r 4 - 8585830400L 14 Mr^ + 29653513376i 12 M 2 r^ - 19705027200L 12 Mr^ 

+ 43981240344 J L 10 M 2 rg - 2392254120(XL 10 Afr[J + 32635169200 J L 8 M 2 r 1 ,° - 16162950800L 8 Mr^ 

+ 12004692860L 6 M 2 r^ 2 - 5726827800L 6 Mr 1 , 3 + 1578230992L 4 Af 2 r 44 - 803087700L 4 Mr^ 5 + 7875L 4 r 46 

- 113069964L 2 M 2 r 1 , 6 + 24092400L 2 Mr^ 7 - 1118250 J L 2 r ( ^ 8 - 13324056M 2 r 1 , 8 + 1526700Mr 1 , 9 + 553875r 20 ) 

+ 2r 2 (L 2 + r 2 ) 3 (l443840000L 18 M 3 - 645120000 J L 18 M 2 r + 3481866240L 16 M 3 r 2 + 376780800L 16 M 2 r^ 

- 857472000L 16 Mr 4 - 5698068736L 14 M 3 r 4 + 12906055680L 14 M 2 r^ - 4598832000L 14 Mr^ 

- 30679784768L 12 M 3 rj5 + 36702979536L 12 M 2 r^ - 10214544200L 12 Mr^ - 46687160912L 10 M 3 rg 
+ 47920180732 J L 10 M 2 r^ - 1203425940(XL 10 Mr 1 1 - 34910457500L 8 M 3 r^° + 33450421620L 8 M 2 r 1 1 1 

- 7955786400L 8 Mr 42 - 13231827540L 6 M 3 r 1 , 2 + 12237529680L 6 M 2 r 1 1 3 - 2830778675L 6 Mr 44 

- 2081061396L 4 M 3 r 1 , 4 + 1920152080L 4 M 2 r^ 5 - 448289925L 4 Mr 1 , 6 - 1750i 4 ^ 7 - 4481148L 2 M 3 r 1 1 6 
+ 24248796L 2 M 2 rf - 11278125L 2 Mrf - 8750L 2 r 1 , 9 + 11067088M 3 r 1 1 8 - 6431148M 2 r 1 , 9 

+ 440325Mr 20 + 77000r 21 ), 



560^ 8 [I? + r 2 )^ (r - 2M) ^ + ^ ' 



F m = — ^ . ,,15/2, ( W + W) ' ^ 5 - 39 ) 

where 

r[6] 



FLi = -28000£ 8 (r - L)(L + r ) (ILL 4 - 74r 2 L 2 + llr 4 ) rf 



+ 50E e (L 2 + r 2 ) (19565^ + 6086A^ 8 - 113785L 2 ^ 7 + 1633964L 2 Mr 1 , 6 + 76615 J L 4 rJ 5 + 2559418L 4 Mr 1 , 4 

- 5075L 6 r^ 3 - 19625630 L 6 Mr] 2 - 100527512L 8 Mr 1 , - 206284326L 10 Mr 8 , - 232489984L 12 Mr^ 

- 151293440L 14 Mr 4 - 53575680L 16 Mr 2 - 8028160L 18 M)r^ 

- E 4 (L 2 + r 2 ) 2 (I094625r 22 + 545450Mr 21 - 4202625L 2 r 20 - 16774936M 2 r 20 + 110 1 1000L 2 MrJ 9 
+ 1414875L 4 ^ 8 - 331621052L 2 M 2 r 1 , 8 - 82244700(XL 4 Mr 1 1 7 - 7875L 6 r 1 , 6 + 1429685188L 4 M 2 r 1 l 6 

- 9480504900L 6 Mr 1 , 5 + 19802086804 J L 6 M 2 rJ 4 - 35145688050L 8 Mr^ 3 + 72068652100L 8 M 2 r 1 1 2 

- 68030273700 J L 10 Mr 1 , 1 + 130518064824L 10 M 2 r 40 - 76873222400L 12 Mr 9 + 129714899808i 12 M 2 r 8 

- 51274604800L 14 Mr^ + 65633972928L 14 M 2 r[j - 18785894400L 16 Mr^ + 8353439744L 16 M 2 4 

- 2924544000L 18 Mr 3 1 - 611471360(XL 18 M 2 r 2 - 2015232000L 20 M 2 )r^ 

+ 2£ 2 (L 2 + rlf (242375r 23 + 138500Mr 22 - 41 9 1 25L 2 r 21 - 10992820M 2 r 21 + 20399392M 3 r 20 
+ 6959650L 2 Mrf + 9625i 4 rJ 9 - 81665756L 2 M 2 r 1 , 9 + 138887552L 2 M 3 rJ 8 - 822884550L 4 Mr 1 , 8 

- 875L 6 rJ 7 + 3518901164L 4 M 2 r^ 7 - 3802035608L 4 M 3 r,i 6 - 7067551050i 6 MrJ 6 + 30648837404L 6 M 2 r 1 , 5 

- 33234129320L 6 M 3 r^ 4 - 25352340750i 8 Mr^ 4 + 106617027800 J L 8 M 2 r^ 3 - 111740075320L 8 M 3 r^ 2 

- 4966003620CLL 10 Mr 1 1 2 + 197830827680L 10 M 2 r 41 - 195029907128L 10 M 3 r 40 - 57553509200L 12 Mr 40 
+ 209550665904L 12 M 2 r[J - 183717663248L 12 M 3 rg - 39556896400L 14 Mrg + 121346652064L 14 M 2 r^ 

- 77791192288L 14 M 3 ^ - 14956032000L 16 Mr^ + 28755231744L 16 M 2 ^ + 7146388480L 16 M 3 r 4 

- 2403072000L 18 Mr 4 - 361973760(XL 18 M 2 r^ + 18731642880L 18 M 3 r 2 - 2297856000 J L 20 M 2 r 
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+ 4902912000L 20 M 3 )r% 

+ (2M - r ) (L 2 + r 2 ) 4 (53375r 23 + 173 6 00Mr 22 + 29750L 2 r 21 - 5212944M 2 r 21 + 10317184M 3 r 20 

- 14183750L 2 Mr 20 + 25375L 4 r 1 l 9 + 62499436L 2 M 2 r 19 - 76153328L 2 M 3 r 18 - 679641900L 4 Mr 18 

+ 7000L 6 r 1 1 7 + 3345821696L 4 M 2 r 17 - 4111469784L 4 Af 3 r 1 , 6 - 5518261350L 6 Mr 46 + 25823273900 J L 6 M 2 r 15 

- 30032966752L 6 M 3 r 14 - 20150264400L 8 Mr 14 + 88797134680L 8 M 2 r 1 l 3 - 96877958296L 8 M 3 r%* 

- 40738068000 J L 10 Mr 12 + 166165277616 J L 10 M 2 r 1 , 1 - 165473393280L 10 Af 3 r 1 l ° - 48857572400L 12 Afr 1 l ° 
+ 177451524416L 12 Af 2 r 9 1 - 149256330784L 12 M 3 rg - 34733020000L 14 Mrg + 101711045376L 14 M 2 r 7 ) 

- 51263318208L 14 M 3 ^ - 13563648000L 16 Mrg + 21232814080L 16 M 2 r^ + 21391316992i 16 M 3 r 4 

- 224716800(XL 18 Mr 4 - 5520998400L 18 M 2 r 3 + 23777402880L 18 M 3 r 2 - 2580480000L 20 M 2 r 
+ 5775360000L 20 M 3 ), 



F% 6] = 875E 8 (-105L 6 + 1189r 2 L 4 - 1531r 4 L 2 + 247^) rf 

- 25E 6 (L 2 + r 2 ) (27055r 47 + 25612Afr 16 - 123235L 2 r 1 1 5 + 1887182L 2 Afr 14 + 62125L 4 r 1 1 3 + 676066L 4 Afr 1 l 2 

- 2625L 6 r 1 , 1 - 26099670L 6 Mr 1 l ° - 90366550L 8 Mrt - 137437024L 10 Mrg - 110937600L 12 Mr 4 

- 46551040L 14 Mr 2 - 8028160L 16 M)r 1 1 3 

+ 2E 4 (L 2 + rlf (370125r 20 + 697975Mr^ 9 - 1065750L 2 r 1 1 8 - 6904028M 2 r 18 + 27977700 L 2 Mr^ 7 

+ 244125L 4 r 1 l 6 - 86371482L 2 M 2 r 1 , 6 - 315079050L 4 Mr 15 + 615225146L 4 M 2 r 14 - 2582807100L 6 Mr 13 
+ 5441132830 J L 6 M 2 r 12 - 7522573725L 8 MrJ 1 + 15199781250L 8 M 2 r 1 l ° - 11253096600L 10 Mrg 
+ 20574272172L 10 M 2 rg - 9303284800L 12 Mr 7 1 + 13728299088L 12 M 2 r[; - 4056729600L 14 Mr^ 
+ 3016609536L 14 Af 2 r 4 - 731136000L 16 Afr 3 1 - 1087846400L 16 M 2 r 2 - 503808000L 18 M 2 )rg 

- E 2 (L 2 + r 2 ) 3 (314125r 21 + 970000Mr 20 - 358750L 2 r 19 - 18222440M 2 r 19 + 31216064APV 1 , 8 

- 3780450L 2 Afr 18 - 875L 4 r 17 - 36108732L 2 M 2 r 1 1 7 + 87877152L 2 M 3 r 46 - 1092352500L 4 Mr 46 
+ 4756636984L 4 M 2 r 15 - 5211718840X 4 M 3 r 44 - 7494234450L 6 Mr 44 + 32418083300 J L 6 M 2 r 13 

- 35018994560L 6 M 3 r 12 - 21674793600L 8 Mr 1 l 2 + 89821855320L 8 M 2 r 1 1 - 92610055960L 8 M 3 r 1 ,° 



- 33236721600L 10 Mr 1 l ° + 127823271040L 10 M 2 r 9 ) - 120667329776L 10 M 3 rg - 28422864400L 12 Mr, 
+ 95019791488£ 12 M 2 r£ - 72612130368L 12 M 3 r[; - 12853344000L 14 Mr^ + 30055494144L 14 M 2 ^ 

- 5260183040 J L 14 M 3 r 4 - 2403072000L 16 Mr 4 - 1609113600L 16 M 2 r 3 + 14441594880L 16 M 3 r 2 

- 2297856000 J L 18 M 2 r + 4902912000L 18 Af 3 )r^ 

+ (r - 2M) (L 2 + r 2 ) 4 (27125r 21 + 299600Afr 20 + 96 25L 2 r 19 - 4164624Af 2 r 1 , 9 + 7521664M 3 r 18 

- 12312300L 2 Afr 18 + 3500L 4 ^ 7 + 60136468L 2 M 2 r 1 , 7 - 77649520L 2 Af 3 r 16 - 435265950L 4 Afr 16 
+ 2136130980L 4 M V 1 . 5 - 2610250432L 4 Af 3 r 1 , 4 - 2918350050L 6 Mr 1 , 4 + 13485739400L 6 M 2 r 13 

- 15471992072L 6 M 3 r 12 - 8700997200L 8 Mr 12 + 37444962000L 8 M 2 r 1 1 1 - 39767055768L 8 M 3 r 1 l ° 

- 13875397200L 10 Alr 1 l ° + 54025898376L 10 Af 2 r^ - 50546368608L 10 M 3 r^ - 12345326000L 12 Mr 8 1 
+ 40319920928L 12 M 2 r 7 - 27561410368L 12 M 3 ^ - 5798688000L 14 Mr|> + 11983523840L 14 M 2 r^ 
+ 2639284736L 14 M 3 r 4 - 1123584000L 16 Alr 4 - 1631539200L 16 Af 2 r 3 1 + 9361981440L 16 M 3 r 2 

- 1290240000L 18 M 2 ro + 2887680000i 18 M 3 )r 2 , 



F m = 0, (5.40) 



F *l = m T ^77t 2^/2 { F *m £ + F mK) , (5.41) 



20 



..(> 

.15 



where 

F$ [6] = 875E 6 rl 3 (l773L 4 r 2 - 337L 6 - 947L 2 r 4 + 15r§) 

+ 175£ 4 r 8 (L 2 + r 2 ) (983040L 18 Af + 7045120L 16 Mr 2 + 21570560L 14 Afr 4 + 36582912L 12 A/r^ 

+ 37139572L 10 Afr^ + 22617566L 8 Mr 1 , + 7708830L 6 MrJ 2 + 225£ 6 r 1 , 3 + 1199730L 4 Afr 1 , 4 - 9375L 4 r 1 , £ 
+ 4926L 2 Air 46 + 9375L 2 r 47 - 225r 49 ) 

+ E 2 rl (L 2 + rlf (2015232000L 20 M 2 + 8400691200L 18 Af 2 r 2 , + 1376256000L 18 Air 3 l + 11063078912L 16 M 2 r 4 
+ 7276953600L 16 Mrg + 82626752L 14 M 2 r^ + 15631481600L 14 M r 7 - 12481032128L 12 M 2 r 8 
+ 17141443200L 12 Air^ - 10632497080£ 10 Af 2 r 1 1 ° + gemieSOOL^Air 1 , 1 - 2120762600L 8 Af 2 r 1 1 2 
+ 2046335900L 8 Afr 1 l 3 + 1035946292L 6 A/ 2 r 1 l 4 - 316454600L 6 Mr 1 , 5 + 15750L 6 ^ 6 + 431941952L 4 Af 2 r 1 , 6 

- 166655300L 4 AfrJ 7 + 133875L 4 ^ 8 + 17346492L 2 A/ 2 r 1 , 8 - 2290400 J L 2 Air 1 , 9 - 850500L 2 ^ 

- 312480Af 2 r 20 + 39375r 22 ) 

- (L 2 + rlf (2580480000L 20 Af 2 r - 5775360000L 20 A/ 3 - 21381242880i 18 Af 3 rg + 4448665600L 18 A/ 2 r 3 , 
+ 2247168000L 18 A/r 4 - 16436058112L 16 APV 4 - 20228515840L 16 M 2 r^ + 12144384000i 16 Mr^ 

+ 37967372736L 14 A/ 3 r^ - 78865174016 J L 14 A/ 2 r 7 l + 27288335200L 14 Mrg + 91117248928L 12 A/ 3 r 8 l 

- 114866020480L 12 A/ 2 r 9 1 + 32738062000L 12 Afr^° + 80974789248L 10 Af 3 r 1 l ° - 86565871136 J L 10 Af 2 r 1 Q 1 
+ 22304895200L 10 Afr 1 , 2 + 35112838392L 8 Af 3 r 1 1 2 - 34540540744L 8 A/ 2 r 1 l 3 + 8392988800L 8 Afr 1 l 4 

+ 6757023136L 6 Af 3 r 1 , 4 - 6373891596L 6 A/ 2 r 1 , 5 + 1516716950L 6 Mr 1 , 6 - 7000L 6 r 1 l 7 + 273916728L 4 M 3 r 1 , 6 

- 286552800L 4 M 2 r 1 1 7 + 80573500L 4 AIr 1 , 8 - 28875LV 1 , 9 - 24711184 J L 2 A/ 3 r 1 l 8 + 14372004L 2 M 2 r 1 , 9 

- SSSOSOZ^Afr 2 , - 50750L 2 r^ + 309120Af 3 r 20 - 24 5 2 80M 2 r 2 , 1 + 131257- 23 ), 

F£ [6] = -875£ 6 r 23 (-105L 6 + 829L 4 r 2 - 587i 2 r 4 + 15rg) 

+ YI5E A rl° (L 2 + r 2 ) ( - 491520L 16 Af - 3092480L 14 Afr 2 - 8102400L 12 Mr 4 - 11336736L 10 Af ^ 

- 8972282 J L 8 A/r 8 ) - 3855372L 6 Mr 1 , - 750282 J L 4 Mr 1 , 2 + 3825LV 1 , 3 - 12696Z 2 Mr 1 , 4 

- 5550L 2 r^ 5 + 225r^ 7 ) 

- E 2 rl (L 2 + r 2 ) 2 (l007616000L 18 Af 2 + 3318681600L 16 Af V 2 , + 688128000L 16 Mrg + 2674925056L 14 M 2 r 4 
+ 3036364800L 14 Afrg - 2164346848L 12 A/ 2 r^ + 5191177600L 12 Afr^ - 4277576360L 10 Af 2 r 8 ) 

+ 4156658800L 10 A/r 9 l - 1698491920L 8 Af 2 r 1 1 ° + 1358613200L 8 Afr^ + 281096500L 6 Af 2 r 1 , 2 

- 46924500L 6 Mrj 3 + 248366216L 4 Af 2 r^ 4 - 97522600L 4 Afr 1 1 5 + 7875LV 1 , 6 + 15819372L 2 Af 2 r 1 1 6 

- 3686900L 2 Mr 1 , 7 - 456750L 2 r 1 l 8 - 312480M V 1 , 8 + 39375r 20 ) 

+ r 2 (L 2 + r 2 ) 3 ( - 2887680000L 18 Af 3 + 1290240000£ 18 Af 2 r - 8163901440L 16 A/ 3 r 2 + 1095372800L 16 M 2 r 3 
+ 1123584000L 16 A/r 4 - 1209975296L 14 A/ 3 r 4 - 11012229120i 14 M 2 r^ + 5089056000L 14 A/r^ 



19718951872L 12 M 3 r£ - 29772000928L 12 M 2 r 7 + 924391 1600L 12 Mr 8 + 28381407744L 10 M 3 r 8 
31906051368L 10 Af 2 r^ + 8496115600L 10 Mr 1 l ° + 16529207256£ 8 M 3 r 1 , - 16531893472L 8 M 2 r 1 , 1 
4060456400L 8 Afr 1 , 2 + 4051974744L 6 Af V 1 , 2 - 382035 1368L 6 M 2 r 1 , 3 + 903318850L 6 Mr^ 4 
233659760L 4 Af 3 r 1 , 4 - 228732252L 4 M 2 r 1 1 5 + 59520650L 4 Air^ 6 - 3500LV 1 , 7 - 17332624L 2 Af 3 r 1 , 6 
10640724L 2 A/ 2 r 1 l 7 - 891800L 2 Mr^ 8 - 11375L 2 r£ 9 + 309120Af 3 r^ 8 - 245280A/ 2 r 1 l 9 + 13125r^). 



D. Electromagnetic case 

In the electromagnetic case, an ambiguity arises in the definition of u a in the angular directions away from the 



world- line. In (2.19) one is free to define u a {x) as they wish provided lim a; _ >2 u a (x) — u a . A natural covariant choice 
would be to define this through parallel transport, u a {x) = g a iu b . However, in reality it is more practical in numerical 
calculations to define u a such that its components in Schwarzschild coordinates are equal to the components of u a in 
Schwarzschild coordinates. Doing so, the regularization parameters are given by 

_ rp sgn(Ar) Er Q sgn(Ar) 

F *[- 1 ]-"2(L2 +r 2)' ^[- 1 ]-2(r -2A/)(L2 + r 2 )' F °™ ~ °' ^ " °' {5A2) 
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F t[0] = f ° /2 (rlJC + 2L 2 £) , (5.43) 

Trro (rg + L 2 ) 

^1 = 3( 2+r2\W ( ^ + ^ ' (5 - 44) 

7rrg (rg + L 2 ) (r - 2M) v ' 



F r £ [0] = 2E 2 L 2 r% + (l? + rl) (2L 2 + r 2 ) (2M - r ), 
^oi=^ + rg (i 2 +^)(r -2M), 



fyo] - 0, (5.45) 



F *[°] = r /TiTT^ [ £ ( 2 ^ 2 + ~ /Cr °] ' (5 - 46) 
TTLr ^L 2 + rf t 

F m = - 2 . 7/2 + ^ - ( 5 - 47 ) 

27rr^ (L 2 + rg) v 7 

2£; 2 ^ (-L 4 + lQL 2 rl + 3r 4 ) + (L 2 + r 2 a ) (60L 6 M + 168L 4 Mr^ + 182L 2 Mr 4 - 13i 2 r^ + 58Mr§ - 5r 7 ) , 
2r 2 (h 2 + r 2 ) (-21L 4 M - 48L 2 Mr 2 + 3L 2 r% - 23Mr 4 + r 5 ) - E 2 r 7 (ILL 2 + 3r 2 ) , 



^ = - 2^(^ + r g^, r „-2M, (^ + ^)' ,5 - 18) 

-2E 4 r 8 (-L 4 + 10L 2 r 2 + 3r 4 ) - (L 2 + r 2 f (2M - r ) (44£ 6 M + 94£ 4 Mr 2 + SAL 2 Mr 4 + L 2 r\ + 3r 7 ) 
+ 2E 2 r% (L 2 + rl) (-30L 6 M - 86i 4 Mr 2 + L 4 rl - 98L 2 Mr 4 + lOL 2 r% - 26Mr§ + r 7 ) , 
E 4 rl Q (Hi 2 + 3r 2 ) - rl (I 2 + r 2 a f (r - 2M) (22L 4 M + 24L 2 Mr 2 , + 2i 2 r^ + 3rfj) 
+ E 2 4 (L 2 + rg) (42L 4 M + 98L 2 Mr 2 - lL 2 r% + 40Mr 4 + r%) , 

F 9[2] - 0, (5.49) 



^-- ^r ^^)' <5 ' 50) 



F 0[ 2 ] = (-2^ 4 - 7i 2 r 2 + 3r 4 ) - (L 2 + r 2 ) (UL 6 M + ML 4 Mr 2 + 54L 2 Mr 4 + L 2 r% + 3r 7 ) , 
Ff[2] = r 2 (L 2 + r 2 ) (22L 4 M + 2AL 2 Mrl + 2L 2 r\ + 3r 5 ) - E 2 r 7 (3r 2 - L 2 ) , 
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where 



Ff w = ~30E 4 r^ (3i 6 - 102L 4 r 2 + 43L 2 r 4 + 20rg) 

+ 2E 2 4 (L 2 + r 2 .) (34560i 12 M + 169728 J L 10 Mr 2 + 333564L 8 Mr 4 + 328912L 6 Mr^ + 167074L 4 Mrg 

- 1245Z, 4 r|] + 32948L 2 Mr 1 ) ° + 1230-Z/VJ 1 + 230M?- 1 , 2 + 555/- 1 , 3 ) 

- 4 (L 2 + rl) 2 (11520L 14 M 2 - 18240L 12 M 2 r 2 + 33600L 12 Mr 3 - 226624L 10 M 2 r 4 + 153960L 10 Mr[> 

- 496164L 8 M 2 r^ + 280764L 8 Mr 7 - 485652L 6 M 2 rg + 255197L 6 Mrl - 230930L 4 M 2 r^° + \\<o77\L 4 Mr\} 

- 30L 4 r^ 2 - 4520(XL 2 M 2 r 1 1 2 + 21487L 2 Mr 1 ? + 270L 2 r 1 1 4 - 1342M 2 r 1 , 4 + 229Mrf + 120r^ 6 ), 
F t f 4] = 15SV 1 , 6 (-87L 4 + 66L 2 r 2 + 25r 4 ) 

- AE 2 r 7 (L 2 + rl) (8640L 10 M + 34872i 8 Mr 2 + 53283L 6 Mr 4 + 37555L 4 Mr^ - 225L 4 r£ + 9809L 2 Mrg 
+ 420L 2 rfJ + 265MrJ° + 165r^) 

+ r 2 (i 2 + rlf (23040L 12 M 2 - 56640L 10 M 2 r 2 + 67200L 10 Mr 3 - 402608L 8 Af 2 r 4 + 249120L 8 Mrl 

- 643736L 6 M 2 r^ + 346608L 6 Mr 7 - 427796L 4 M 2 rg + 217192L 4 M^ - 110916L 2 M 2 ^ 
+ 52580L 2 Mr 1 , 1 + 615L 2 r^ 2 - 5368M 2 r^ 2 + 1516Mr^ 3 + 2557- 1 , 4 ), 



if [4] = 30£ 6 ^ 7 (3L 6 - 102i 4 r 2 + 43i 2 r 4 + 20rg) 

- E 4 rl (L 2 + r 2 ) (69120L 12 M + 339456L 10 Mr 2 + 667128L 8 Mr 4 + 657884L 6 Mrl ~ SOL 6 r 7 + 334658£ 4 Mrg 

- 2745L 4 r[J + 62656L 2 Mrj" + 4080 J L 2 r* 1 + 610Mr^ 2 + 1035r^ 3 ) 

+ (i 2 + rl) 3 (2M - r )(46080i 14 M 2 + 89856L 12 M 2 r 2 + 53760L 12 Mr 3 1 - 86336L 10 M 2 r 4 + 211520L 10 Mr^ 

- 344128L 8 M 2 r^ + 317600L 8 Mr 7 - 306808L 6 M 2 r^ + 221140L 6 Mr^ - 98676L 4 M 2 r^° + 66220£ 4 Mr 1 , 1 
+ 60LV 1 , 2 - 5244 J L 2 M 2 r^ 2 + 4440L 2 Mr 1 , 3 + 105L 2 r^ 4 + 160M 2 r^ 4 - 75r^ 6 ) 

+ 2E 2 r\ (L 2 + rl) 2 (23040L 14 M 2 - 36480L 12 M 2 r 2 + 67200L 12 Mr% - 445056L 10 M 2 r 4 + 303824i 10 Mr^ 

- 959952L 8 M 2 rg + 545340L 8 Mr^ - 922996L 6 M 2 rg + 486240L 6 Mrl - 428644 J L 4 M 2 r 1 1 ° + 216604 J L 4 Mr 1 1 1 
+ 105L 4 ^ 2 - 78428L 2 M 2 r 1 , 2 + 35368L 2 MrJ 3 + 1350L 2 r^ 4 - 2684M 2 r 1 , 4 + 608Mr^ 5 + 165r^ 6 ), 

F^ [4] = -lbE G r™ (-87L 4 + 66L 2 r 2 + 25r 4 ) 

+ E 4 r^ (L 2 + rl) (34560L 10 M + 139488L 8 Mr 2 + 213132L 6 Mr 4 + 150250L 4 Mr£ - 9lbL 4 r 7 Q + 37496£ 2 Mrg 
+ 2bhOL 2 rl + 1210Mr^° + 585T- 1 , 1 ) 

- rl (L 2 + rl) 3 (r - 2M) (63760L 8 M 2 r 4 - 2304(XL 12 M 2 - 24768L 10 M 2 r 2 - 26880L 10 Mr 3 ) - 82240L 8 Mrl 
+ 115848L 6 M 2 rg - 87920L 6 Mr 7 l + 56084L 4 M 2 rg - 36340L 4 Mrl + 5324L 2 Af 2 r 1 ) ° - 3540L 2 Mr 1 , 1 

+ 15i 2 ri 2 - 160M 2 ri 2 + 75^) 

- E 2 rl (L 2 + rlf (23040L 12 M 2 - 56640L 10 M 2 r 2 + 67200L 10 Mr 3 - 394416L 8 Af 2 r 4 + 245024L 8 Mrl 

- 618528L 6 M 2 rg + 334004L 6 Mr 7 - 400636L 4 M 2 r 8 + 203252L 4 Mr^ + lSOL 4 ^ - 97892L 2 M 2 r 1 [ ° 
+ 44568i 2 Mr^ + l365L 2 rl 2 - 5368M 2 r^ 2 + 1816MrJ 3 + 105r^ 4 ), 




(5.52) 



where 



^[4] = 0, 



(5.53) 



F m = 



AQuLrl 1 {IP- + rl) 



(5.54) 
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where 

F| [4] = -VoE^rl 4 (2L 6 + \lh\ 2 - 108L 2 r 4 + 5r^) 

+ 2E 2 r 7 (L 2 + r 2 ) (4096L 10 Af + 16188L 8 Air 2 + 24154L 6 Afr 4 + 16608L 4 Afr^ - l65L 4 r 7 + 5986L 2 Mr% 

- 810i 2 rg + 75r^) 

+ (L 2 + rlf (46080L 14 A/ 2 + 89856L 12 Af 2 r 2 + 5376(LL 12 A/rg - 86336L 10 Af 2 r 4 + 211520L 10 A/rg 

- 344128L 8 A/ 2 r^ + 317600L 8 Mr 7 - 306808i 6 Af 2 rg + 221140£ 6 Afr^ - 98676L 4 Af V 1 , 

+ 66220L 4 Afr 1 , 1 + 60L 4 ^ 2 - 5244L 2 A/ 2 r^ 2 + 4440L 2 Mr 1 , 3 + 105L 2 r^ 4 + 160A/ 2 r^ 4 - 75r^ 6 ), 
Ff w = 15FV 1 , 6 (i 4 - 58L 2 r 2 + 5r 4 ) 

- 2F 2 r[J (L 2 + r 2 ) (2048L 8 A/ + 6302L 6 Af^ + 6790L 4 ATr 4 - 90i 4 rg + 3256L 2 Air^ - 375L 2 r£ + 75rg) 
+r 2 (L 2 + r 2 ) 2 ( - 23040L 12 A/ 2 - 24768L 10 Af 2 r 2 - 26880L 10 Afr 3 + 63760L 8 Af 2 r 4 - 82240 J L 8 Afrg 

+ 115848L 6 Af 2 rg - 87920L 6 Afr^ + 56084L 4 M 2 rg - 36340L 4 Afrg + 5324L 2 Af 2 r 1 1 ° - 3540L 2 A:fr 1 1 1 

+ \hL 2 r 12 - VS)M 2 r 12 + 75r 1 , 4 ). 

E. Gravitational case 

1. Self-force regularization 

The self force on a gravitational particle is given by 

pa = k abcd hf c . d , (5.55) 

where 

k abcd = l g ad u b u c _ g ab u c u d hf^^d + I u »^ u «« + ^g ad g bc . (5.56) 

Note that, as in the electromagnetic case, an ambiguity arises here due to the presense of terms involving the four- 
velocity at x. One is free to arbitrarily choose how to define this provided \im x ^ s u a — u a . Following Barack and 
Sago [S3], we choose to take the Schwarzschild components of the four velocity at x to be exactly those at x. The 
regularisation parameters in the gravitational case are given by 

^ 1 l =T 2(L2 + r g°)(° _2M)' F [- 1 ] =T 2(L2+ r 2)' F K]= ' F K1 = °' ^ 



Fmi 



nrK^ + r 2 )^ { ^ ~ + + ^ " ^^O^ " r ° ^ " + r ") + 



'•o 



^[o]=0, F* =- - 7377 ^- ^ 1/9 [(2£ 2 +r 2 )g-r 2 /C], (5.59) 



7rir3(i 2 



[(2I/ + ^)£-r ^C], 



^ = 2^ (n , - 2A^ + ^ ^ + *W 



where 



F* [2] = -2F 2 r 5 (llL 4 + 34L 2 r 2 + 15r 4 ) 

-(L 2 + r 2 )(276L 6 Af + 768L 4 A/r^ + 782L 2 Afr 4 - 37L 2 r^ + 274Afr[j - 29r£), 
F£ [2] = E 2 r%(YlL^ + 35L 2 r 2 + 15r 4 ) + 2r 2 (L 2 + r 2 )(93L 4 M + 204L 2 A/r 2 - 9L 2 r^ + 107Afr 4 - 7r%), 
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2nr 7 (L? + r 2 ) 7 / 

where 



F k = oZJFTWT^mW + F mV, (5-61) 



F £[2] = ( r o - 2M)(L 2 + r 2 ) 2 (188L 6 M + 406L 4 Mr 2 + 222L 2 Mr 4 + 13L 2 r§ + 15r 7 ) - 2£ 4 r 8 (llL 4 + 34L 2 r 2 - 
- 2£ 2 r 3 (L 2 + r 2 )(138L 6 M + 422L 4 Mr 2 - 19L 4 r 3 + 422L 2 Mr 4 - 34L 2 r^ + \22Mr% - 7r 7 ), 

F K[2] = - r o( r o - 2M)(L 2 + r 2 ) 2 (94L 4 M + 96L 2 Mr 2 + 14L 2 rg + 15r 5 ) + £ 4 r 8 (12L 4 + 35L 2 r 2 + 15r 4 ) 
+E 2 rl{L 2 + r 2 )(210L 4 Af - 12L 4 r + 410L 2 Mr 2 - 19L 2 rg + 184Mr 4 + r%)}, 



15r 4 ) 



where 



if 2] = 0, (5.62) 



F P] ~ 2wL^(L 2 + r 2 )5/2 ( F £ [2] f + F K[2f)' ( 5 - 63 ) 



^£[21 = -(L 2 + rl)(lS8L 8 M + 406L 6 Mr 2 - 64L 6 r 3 + 222L 4 Mr 4 - 163L 4 ^ - 145L 2 r^ - 48r{>) 



£ 2 L 2 r^(38L 4 + 31L 2 r 2 - 15r 4 ), 



F* [2] = E 2 L 2 4(12L 4 + L 2 rl - 15r 4 ) + rg(L 2 + r 2 )(94L 6 M + 96L 4 Mr 2 - 74L 4 rg - 121L 2 rg - 48r£), 



F W 40^(r -2M)V 2 +r 2 )iV2^W f + ^W^). ( 5 - 64 ) 

where 

F| [4] = 30£; 4 ^ (64L 10 + 384L 8 r§ + 989L 6 r 4 + Y222L A r% + 437L 2 rg + 12r^°) 

+ 2E 2 r\(L 2 + rl){92mOL 12 M + 445008L 10 Mr 2 + 859044L 8 Mr 4 - 1920L 8 r^ + 838312L 6 M^ - 9780L 6 r 7 

+ 433114L 4 Mrg - 19065L 4 r^ + 102188L 2 Mr 1 , - 9870Z/V 1 , 1 + 5030Mr^ 2 - 585r^ 3 ) 
+ 4(L 2 + r^) 2 (46080L 14 M 2 + 403200L 12 M 2 rg - 48000L 12 Mr% + 1231984L 10 M 2 r 4 - 219840L 10 Mr^ 
+ 1841004L 8 M 2 r[> - 411324L 8 Mr 7 + 1490772L 6 M 2 r 8 - 406397L 6 Mr^ + 480L 6 r™ + 668810X 4 M 2 r 40 
-232331L 4 Mr 41 + 2040L 4 ri 2 + 161840L 2 M 2 r 42 - 75367L 2 Mr 1 , 3 + 1590L 2 r^ 4 + 18382M 2 r 44 
- 10669Mr^ 5 +210r^ 6 ), 
F K.{i] = -15 J B 4 rJ 2 (64L 8 + 328L 6 r 2 + 495L 4 r 4 + 22£ 2 r^ - 81r§) 

- 4E 2 r 7 (L 2 + rg)(2592(XL 10 M + 102372L 8 Mr 2 + 152523L 6 Mr 4 - 480L 6 r^ + 103375L 4 Mrg - 1575L 4 r£ 
+ 25889L 2 Mrg - 120i 2 r^ - 455Mr^° + 495r^) 

- rl{L 2 + r^) 2 (92160i 12 M 2 + 766080L 10 M 2 r 2 - 130560L 10 Mrj] + 2014208 J L 8 M 2 r 4 , - 497760L 8 Mrl 

+ 2443016L 6 M 2 r|> - 743088L 6 Mr 7 + 1527236L 4 M 2 r 8 - 552712L 4 Mr^ + 960L 4 r™ + 496596L 2 M 2 r 1 ) ° 



206180L 2 MrJ 1 - 135L 2 r 42 + 73528M 2 r 1 1 2 - 30796M^ 3 - 735r, 
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where 

F r m = 30£ 6 r^ 3 (64L 10 + 384L 8 r 2 + 989L 6 r 4 + 1222L 4 r^ + 437L 2 r 8 + Ur] ) 

+ E A rl{L 2 + r 2 )(184320L 12 M + 893856L 10 Mr 2 - 1920L 10 r 3 + 1746888L 8 Mr* - 18240i 8 r^ 

+1744604 J L 6 Mr^ - 53550L 6 rJ + 907298L 4 Mr 8 , ~ 58665L 4 ^ + 197536L 2 MrJ° - 16320L 2 rJ 1 + 9010Mr^ 2 
-645r^ 3 ) 

+ 2E 2 rl{L 2 + r^) 2 (92160i 14 M 2 + 1036800L 12 M 2 r^ - 211200L 12 Mr^ + 3363456L 10 M 2 r 4 - 889424L 10 Mr^ 
+ 5003472L 8 M 2 rg - 1487220L 8 Mr£ + 192(XL 8 ^ + 3866356L 6 M 2 r^ - 1272840i 6 M^ + 9780L 6 r^° 
+ 1570564L 4 M 2 ri° - 594724L 4 Mr 1 l 1 + 10875L 4 ^ 2 + 321308L 2 M 2 r 42 - 146848L 2 Mr 1 , 3 + 1830L 2 rJ 4 
+ 36764M 2 r 1 , 4 - 20288Mr 1 l 5 - 105rJ 6 ) 

- (r - 2M)(L 2 + r 2 ,) 3 (18432(XL 14 M 2 + 1711104 J L 12 M 2 r 2 - 245760L 12 Mr 3 1 + 4872896i 10 M 2 r 4 

- 884480L 10 Mr^ + 6311728i 8 M 2 r[j - 1185120Z, 8 Mr£ + 4083688L 6 M 2 rg - 721620L 6 Mrf| + 1920i 6 r^° 
+ 1299396i 4 M 2 r^° - 198700 J L 4 Mr 1 l 1 + 2460£ 4 r 1 , 2 + 209484L 2 M 2 rf - 28120L 2 Mr^ 3 - 105L 2 r^ 4 

+ 28640M 2 r 1 , 4 - 5120Mr 45 - 525r 46 ), 
F r m = -VSE & rf{<aAL & + 328L 6 r 2 + 495i 4 r 4 + 22L 2 r^ - 81r 8 ) 

- E 4 r™(L 2 + r 2 )(103680L 10 M + 415248L 8 Mr 2 - 2880L 8 r 3 + 627612i 6 Mr 4 - 10680L 6 ^ + 418570L 4 Mr^ 

- 8835L 4 r 7 + 93896L 2 Mrl + 4350L 2 r[J - 2870Mr^° + 2505r") 

- E 2 rl{L 2 + r 2 ) 2 (92160£ 12 M 2 + 1019520£ 10 M 2 r 2 l - 25728(XL 10 Mr 3 1 + 2806176L 8 M 2 r 4 - 893744L 8 Mrl 
+ 3309168£ 6 M 2 rg - 1180004i 6 Mr^ + 1920Z, 6 rj] + 1878796L 4 M 2 rg - 724772L 4 MrfJ - 900LV 1 , 

+ 517652L 2 M 2 rl° - 205608L 2 Mr 1 , 1 - 5685L 2 r^ 2 + 73528M 2 r 1 1 2 - 28696Mr 1 , 3 - 1785r£ 4 ) 
+ r 2 (r Q - 2M)(L 2 + r 2 ) 3 (92160L 12 M 2 + 815232L 10 M 2 r 2 - 157440L 10 Mr 3 + 1939840L 8 M 2 r 4 

- 468160L 8 Mr^ + 1942488i 6 M 2 r^ - 494560L 6 Mr\ + 892484L 4 M 2 r 8 - 217780L 4 Mr^ - 960i 4 r*° 
+ 195164L 2 M 2 r^° - 38820L 2 Mr 1 , 1 - 1545L 2 r^ 2 + 28640M 2 r 1 l 2 - 5120MrJ 3 - 525rJ 4 ), 



Ft = 0, (5.66) 



F [A\ - 4 07rr 13 L 5 (L 2 + ,,2)9/2 (^V + F Uf^ ^ 5 - 67 ) 

where 

F* [4] = 15E 4 L 4 rl°(128L 10 + 960i 8 r 2 + 2266i 6 r 4 + 1369L 4 r^ - 228L 2 r 8 - 35rJ°) 

+ 2E 2 L 2 rl(L 2 + r^)(115200L 14 M + 449744L 12 Afr^ + 660732i 10 Afr 4 - 5760L 10 r^ + 442406L 8 Mrl 

- 21300L 8 r 7 + 116472L 6 Mr 8 - 18475L 6 r|] - ll8SL 4 Mrl° + 431(XL 4 r 1 l 1 + 10765L 2 r 43 + 4240^ ) 

+ (L 2 + r 2 ) 2 (184320L 18 M 2 + 1711104L 16 M 2 r 2 - 24576(XL 16 Mr 3 + 4872896L 14 M 2 r 4 - 884480L 14 Mr£ 
+ 6311728L 12 M 2 rg - 1090720L 12 Mrl + 4083688L 10 M 2 rg - 436180i 10 M^ - 5760 J L 10 r^° 
+ 1299396L 8 M 2 r^° + 73140L 8 Mr 1 , 1 - 83700L 8 r 1 l 2 + 209484L 6 M 2 r 1 l 2 + 15720L 6 Mr 1 ^ - 236585L 6 r^ 4 
+ 2864(XL 4 M 2 r 1 1 4 - 63200L 4 M^ 5 - 271325L 4 rJ 6 - 21120i 2 Mr^ 7 - 138400i 2 r^ 8 - 25600r 20 ), 
F* [4] = -15 J B 4 J L 4 r^ 2 (192i 8 + 584L 6 r 2 + 169L 4 r 4 - 322L 2 r^ - 35rg) 

- 2E 2 L 2 r 7 Q (L 2 + r 2 )(63360L 12 M + 197992i 10 Mr 2 + 216538i 8 Mr 4 - 4800i 8 r^ + 87890L 6 Mr^ - 7U0L 6 r 7 Q 
+ 5264i 4 Mrg + 2455£ 4 rfJ + 8645i 2 r^ + 4240r^ 3 ) 

- r 2 (L 2 + r^) 2 (92160L 16 M 2 + 815232£ 14 M 2 r 2 , - 157440i 14 Mr^ + 1939840L 12 M 2 r 4 - 422080i 12 Mr^ 

+ 1942488L 10 M 2 r^ - 309120£ 10 Mr£ + 892484L 8 M 2 r 8 + 9580L 8 Mrl - 39360L 8 r*° + 195164L 6 M 2 r^ 
+ 22780L 6 Mr 41 - 152745 J L 6 r* 2 + 28640L 4 M 2 r 42 - 52640X 4 Mr* 3 - 213325i 4 r* 4 - 21120L 2 MrJ 5 

- 125600L 2 ^ 6 - 25600r 1 , 8 ). 
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2. huu regularization 



The quantity 



(5.68) 



was first proposed by Detweiler [S] as a tool for constructing gauge invariant measurements from self-force calculations. 
It has since proven invaluable in extracting gauge invariant results from gauge dependent self- force calculations [61U77) . 

Much the same as with self-force calculations, the calculation of H^ R ~> requires the subtraction of the appropriate 
singular piece, u a u from the full retarded field. In this section, we give this subtraction in the form of 

mode-sum regularization parameters. In doing so, we keep with our convention that the term proportional to I + | 
is denoted by (= in this case), the constant term is denoted by -ff[o], and so on. 

Note that, as in the self-force case, an ambiguity arises here due to the presense of terms involving the four-velocity 
at x. One is free to arbitrarily choose how to define this provided lim a; _ i . s u a — u a . As before, we choose this in such a 
way that the Schwarzschild components of the four velocity at x are exactly those at x. The regularisation parameters 
are then given by 



2K 



[o] 



(5.69) 



H 



[2] 



r3(L 2 + r 2 )3/ 2 ' 



(5.70) 



where 



HL = 2E 2 r% 



(L 2 

Hf 2] = -£ 2 r 3 (16L 2 



rl){mL 2 M~ 
- 17r§) - 2(L 2 



8L 2 r 



38Mr, 



o 



r 2 )(16L 2 M 



9r 3 ) 
4L 2 r 



33Mr 2 - I2r%), 



20nr w (L 2 + r 2 ) 7 / 2 



(5.71) 



where 



Hfc = -120E 4 



ttK, 

H [4] 



1 {8L 4 



+ 120L 4 r^ - 



YIL 2 rl + 7r 4 ) + 2E 2 r%{L 2 + r 2 )(3584L 8 Af + 



!27l2L 6 Mr 2 



15516L 4 Afr 4 



2(L 2 



f 6182L 2 Mr% + 735L 
(1536L 10 M 
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13888L 8 M 2 r 2 



34Mr 8 l + 4954) 

1600L 8 Mrg + 40584L 6 M 2 r 



9440L 6 Mrjj 



46888i 4 M 2 ? 



15S 4 ri° 



U100L 4 Mr 7 + 120£ 4 r^ + 18936L 2 M 2 r^ 
(64L 6 + 224L 4 r 2 + 259i 2 r 4 + 91rg) 



5350L 2 Mrl + 15L 2 r^° + 340M 2 r 1 , + 850Mr^ - 90 1 , 2 ), 



3174L 4 Mr 2 



- 4E 2 r 7 Q (L 2 + r 2 )(1376L 6 M 

- r 2 {L 2 + r 2 ) 2 (1536i 8 M 2 + 15904L 6 M 2 r 2 



420L 4 r2 



1965L 2 Mr 4 



960L 2 ^ 



221Mr% + 510rJ) 



- 11040L 2 Mr 7 ~ 720L 2 r^ 



680M' 



7360L 6 Mrl 
860Mr2 - 705r, 



36160L 4 Af 2 r 4 - 19320L 4 Mr, 



22412I/M 



2_6 
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F. Example 

As an example application of our high order regularization parameters, we consider the case of a scalar particle on a 
circular geodesic of the Schwarzschild spacetime. In this case, the retarded field may be computed using the frequency 
domain method described in |53j , along with improved asymptotics for the boundary conditions (by expanding inside 
the exponential rather than outside) and with the use of the arbitrary precision differential equation solving support 
in Mathematica |78j . These improvements allowed us to substantially increase the accuracy of the computed retarded 
field. We found this to be necessary to get the full benefit from the the higher order regularization parameters. 

In Fig. [2j we show the effect of subtracting in turn the cumulative sums of the regularization parameters from the 
full retarded field. Shown in the figure in order from top to bottom are F* et and the result of subtracting from it in 
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FIG. 2. Regularization of the radial component of the self-force for the case of a scalar particle on a circular geodesic of radius 
ro = 10A/ in Schwarzschild spacetime. In decreasing slope the above lines represent the unregularised self-force, self-force 
regularised by subtracting from it in turn the cumulative sum of A l r , B l r , D l r , Fl, H l r , Jl, L l r , N l r . 



turn the cumulative sum of the regularization terms A l r , B l r , D l r , F^, H l r , J', L l r and N l r . The parameters A r , B r , 
D r , F r and H r are the analytically derived ones given in Sec. V C| The parameters J r , L r and H r were determined 
through a numerical fit to the data. The resulting rapid convergence with I enables the calculation of an extremely 
accurate value for the self-force. Summing over I, we find F r = 0.000013784482575667959(3), where the uncertainty 
in the last digit is estimated by assuming that the error comes purely from the fact that the sum is only done up to 
a finite l max = 80. 



VI. EFFECTIVE SOURCE 



As another application of our high order expansions of the Detweiler- Whiting singular field, we consider its use in 
the effective source approach to calculating the self-force. The effective source approach - independently proposed 
by Barack and Golbourn [37] and by Vega and Detweiler - relies on knowledge of the singular field to derive 
an equation for a regularized field that gives the self-force without any need for post-processed regularization. If the 
singular field is known exactly, then the regularized field is totally regular and is a solution of the homogeneous wave 
equation. In reality, exact expressions for the singular field can only be obtained for very simple spacetimes. More 
generally, the best one can do is an approximation such as that given in Sec. 
approximate singular and regularized parts, 



III Splitting the retarded field into 



(6.1) 



substituting into the wave equation, Eq. (2.1 ), and rearranging, we obtain an equation for the regularized field, 

T>A ~B 

with an effective source, 



S. 



effi 



•2) 



cA 
deft 



s^fs) -4ttQ / u A 5 4 {x,z{T r ))dT'. 



(6.3) 



For sufficiently good approximations to the singular field, (pf 



and S J 



ama u are finite everywhere, in particular on the 
world-line. As a result, one never encounters problematic singularities or (5-functions, making the approach particularly 
suitable for use in time domain numerical simulations. A detailed review of this approach can be found in |79[ I80j . 

In Figs. |3|4[ [5] and [6] we show the result of applying our expansions to the case of a scalar particle on a circular 
geodesic of radius ro = 10M in Schwarzschild spacetime. Similar plots can be obtained for the electromagnetic case, 
gravitational case and for more generic motion. However, the general structure does not change and is best illustrated 
by this simple example. 
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FIG. 3. Terms in the coordinate expansion of the singular field for 0(e 1 ) (top left) to 0(e 6 ) (bottom right). Shown is the 
case of a circular geodesic of radius ro = 10M in Schwarzschild spacetime. 




FIG. 4. Terms in the coordinate expansion of the singular field for C(e *) (top left) to 0(e 6 ) (bottom right). Shown is the 
case of a circular geodesic of radius ro = 10M in Schwarzschild spacetime. 



VII. DISCUSSION 

In this paper we have developed high order expansions of the Detweiler- Whiting singular field of a point scalar 
or electromagnetic particle and of a point mass. Many of our expressions are very general, not necessarily being 
restricted to any particular choice of spacetime. In our explicit coordinate calculations, however, for simplicity we 




FIG. 5. Effective source for the approximate singular field of order 0(e x ) (top left) to C(e 6 ) (bottom right). Shown is the 
case of a circular geodesic of radius ro = 10M in Schwarzschild spacetime. 




chose to limit ourselves to the Schwarzschild geometry. A logical extension of this work would be to consider, instead, 
the Kerr spacetime. This has already been done to 0(e 2 ) for the case of a scalar charge in a circular geodesic orbit in 
Kerr spacetime and spherical- harmonic regularization parameters have been computed [81) . We will present a more 
general version of this calculation in forthcoming work [82] . However, the use of spherical harmonics is not well suited 
to the Kerr spacetime. A more appropriate choice of basis functions are the spheroidal harmonics; it may be more 
sensible to compute regularization parameters in this spheroidal harmonic basis. Nevertheless, the full singular field 
is the same, independently of the choice of basis functions, and the majority of our calculations can easily be carried 
over to the Kerr case with little modification. 

In the cases of electromagnetic and gravitational fields, there is an arbitraryness in the choice of gauge. For 
this work, we have focused only on the case of Lorenz gauge, A a [a — and h a t, ]b = 0, where the singular field 
is best understood. However, other gauges such as Regge- Wheeler and Radiation gauge are better suited to time 
domain numerical calculations. Indeed, there are many outstanding problems in the time domain evolution of metric 
perturbations in Lorenz gauge [53]. Given recent developments in understanding the singular behavior of the retarded 
field in other gauges [721 IZ21 Ell [85], it may soon be possible to implement calculations similar to what we have 
presented here for other gauges. This could be achieved, for example, by either taking our Lorenz gauge expressions 
and performing the gauge transformation to a different gauge or by directly computing the singular field in the other 
gauge. It should be noted, however, that there are many non-trivial issues which can arise when considering other 
gauges; for example, the gauge tranformation itself may be singular or the singular field may be even more singula^ 
than in Lorenz gauge. Any calculations will require considerable care to correctly deal with such issues. 

In Sec. VI we gave plots of an effective source which may be used in a numerical evolution. We have not yet, 
however, fed this effective source into an actual numerical code. Existing results from effective source calculations 
have used a source obtained by taking the first four terms in the expansion of the singular field (to C(e 2 )) which yields 
a source which is C . This is sufficient to obtain correct results, but limits the convergence of any numerical scheme 
with resolution or number of modes required. By including the next four terms in the expansion of the singular field 



(to C(e 6 )) the highest order effective source shown in Sec. VI is now C 4 , allowing for the use of high-order finite 
differencing to potentially obtain a significant increase in numerical accuracy. 

We have limited the focus of the present work to linear order. However, for calculations with sufficient accuracy to 
match the requirements of gravitational wave detection, it will be necessary to go beyond first order and to compute 
a second order self-force. Rapid progress is being made in resolving the many conceptual issues which arise when 
trying to derive the second order equations of motion [2"9Tf5T] and the computation of a second order effective source. 
The next step is to implement this second order prescription in a numerical calculation. As this involves an effective 
source which contains terms involving the first order singular field, it seems likely that many of the calculations done 
here may be of use for second order. 



2 A specific example of this is the case of metric perturbations in Regge- Wheeler gauge, where there are additional delta-function 
divergences in the singular field which are not smoothed out by a mode decomposition |71| . 



30 



ACKNOWLEDGEMENTS 

We are extremely grateful to Niels Warburton, Sarp Akcay and Leor Barack for making available their data for 
the spherical harmonic modes of the retarded field in scalar and gravitational cases, and to Roland Hass for making 
the equivalent data available for the electromagnetic case. We thank Sam Dolan, Marc Casals, Jose Luis Jaramillo, 
Michael Jasiulek, Abraham Harte and particularly Ian Vega for many insightful discussions during the progress of 
this work. Finally, we thank participants of the 2010 and 2011 Capra meetings (in Waterloo and Southampton, 
respectively) for many illuminating conversations. 



Appendix A: Coordinate Expansions in Schwarschild Space-Time 



In this Appendix we give coordinate expansions of the key quantities appearing in the singular field, Eqs. (2.131, 



(2.18) and (2.23) 



Synge World function 

Letting cos 7 = cos 9 cos 8' — smOsinO' cos(4> — 0') so that 2(1 — cos 7) = Aw 2 + Aw 2, the expansion of the world 
function to the order required in this paper is 

i+j+2k<9 

a(x,x')= a ^(t' -t)V -rf (l-cos 7 ) fe +0(e 10 ), (Al) 

i,j,/c— 

where the non-zero coefficients are 

, Mr 1 M(UM 2 - 19Afr + 30r 2 ) M 1 

0001 = r , CT002 = 0003 = 7^ M (9r - 2M), a 00i = — , 0020 = ~, 

3 90 840?' r 2 

M(r-2M) M (lOM 2 — HMr + 3r 2 ) M (-92A/ 3 + 142M 2 r - 78Afr 2 + 15r 3 ) 
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6r 2 » — 60r 3 ' — 840r 4 

M 2 {2M - r) M 2 (54A/ 2 - 49Mr + llr 2 ) 



0040 — 7T~. E ! 0041 



24r 5 ' u ^ 360r 6 



M 2 (1956M 3 - 2702M 2 r + 1228Mr 2 - 183r 3 ) M 3 (26M 2 - 25Mr + 6r 2 ) 



0042 — inner, 7 ) 0060 



10080r 7 ' UDU 720r 9 

M 3 (-1818AT 3 + 2475A/ 2 r - 1115Air 2 + 166r 3 ) Af 4 (978A/ 3 - 1393M 2 r + 660A/> 2 - 104r 



3\ 



0061 — : r , f -» ln ) 0080 



15120r 10 ' uou 40320r 13 

M M 1 / ic 7A/ 2 \ M M(4M-r) 

0101 = r, CT102 = -pr, 0103 = 0104 = M 15 — , 0120 = -tt^> 0121 — 



g , -iiM 20 , "xu« r2 y, "izu 2r2 , 12r3 , 

Af (30Af 2 - 22Air + 3r 2 ) M (368Af 3 - 426A/ 2 r + 156A/r 2 - 15r 3 ) Af 2 (2r - 5M) 

CTl22= 120^ ' CTl23 = 1680^ ' aU0= 24^ ■ 
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M 3 (23931M 3 - 31560M 2 r + 13652Mr 2 - 1930r 3 ) M M 
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Van Vleck determinant 



Inserting the above expansion for a(x,x') into the definition of the Van Vleck-Morette determinant, Eq. (3.2), gives 

i+j+2k<7 

A 1 / 2 (x,x') = l+ J2 A l/fe(*' -ty(r'~ry(l- cos 7 ) fc + 0( e 10 ), (A3) 

i,j : k— 

where the non-zero coefficients are 
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400 60r 4 (r-2M) 2 ' 401 50 40r 5 (r - 2M) 3 

1/2 _ Af 2 (-3526M 2 + 3746M r - 981r 2 ) 1/a _ Af 2 (4Af - 3r) 



12,1 5040r 8 (r - 2A/) 2 ' 500 60r 5 (r - 2Af) 3 ' 

^501 — 



i/ 2 M 2 (-1190Af 3 + 2664M 2 r - 2035Afr 2 + 537r 3N > 



5040r 6 (r - 2M) 4 

1/2 M 2 (9392A/ 3 - 15628A/ 2 r + 8631Mr 2 - 1572r 3 ) 1/2 Af 2 (5432Af 2 - 7720A4V + 2943r 2 ^ 



VJ!I 5040r 9 (r - 2M) 3 ' 600 30240r 6 (r - 2A/) 4 

, /2 AT 2 ( 1036A/ 3 - 2120M 2 r + 1524Afr 2 - 393r 3 ) 

A 700 = " ^ A4) 



Expansions of an arbitrary point on the world line about x a 

The four velocity of a general geodesic orbit taken to lie in the equatorial plane is given by the standard expres- 
sions [86] 



It is straightforward to calculate the higher order proper time derivatives of these expressions and evaluate both the 
four velocity and its higher derivatives at x a , giving, for example, 

Er /Z ( 2AA ( ■ ■ L 

r -2M y V r o J V r o J r o 

, 2EMr .. L 2 r -Mr 2 -3L 2 M ■■ , 2Lf Q 

to = ~T r 77,72' r °- 74 ' ^ = 73^' 



33 



to 



2EM [2 (E 2 -l)r%- r 2 (3L 2 + 2M 2 ) + 9L 2 Mr + 5Mrg - 6L 2 M 2 ] 

4 (ro - 2Mf 
f (-3L 2 r + 2Afr 2 + Y2L 2 M) 



ro = 



0o = O, 0c 



2L [3 (£ 2 - l) rg - 4L 2 r + IMrl + 9L 2 M] 



(A6) 



Combining Eq. (A6) with Eq. (3.17), we can express x a in terms of x a and At: 



Em 



:At- 



EMr, 



o 



r - 2M ( ro - 2M)' 



rAr 2 + • • • , r' = ro + tqAt — 



-L 2 r + A/r 2 + 3L 2 M) 



2r\ 



An 



7T 

2' 



L ■ Lr A ,2 



rAr 



-Ar^ 



(A7) 



It is also straightforward to obtain 8x a , in terms of Ax a , x a and At by noting that Sx a — x a — Ax a — x a . Finally, 
we can calculate u a in terms of x a and At by inserting r' from Eq. ( A7 ) into our equations for the four velocity: 



Er 2EMrf 



At 



r - 2M ( ro _ 2M)' 

EM {6L 2 M 2 - 9L 2 Mr + rl [3L 2 + 2M 2 - 2 (E 2 - l) rl - 5Mr ] } 
+ 4 (2M - r ) 3 



At 2 



u = r + 



r Q (L 2 - Afr ) - 3i 2 M r (l2L 2 M - 3L 2 r + 2Afr 2 ) 

-At 1 



2rg 



At 2 



if = 0, 



, _ L 2Lr , L [3 (E 2 - l) rg - 4L 2 r + 7Afr 2 + 9L 2 A/] 



r 2 r^ 



A7 



(A8) 



Expansions of retarded and advanced points 

Taking At to have leading order e, the same leading order of our Ax terms, we can further expand it in orders of 
e. giving 



At = Tie + T 2 e 2 + T 3 e 3 + T 4 e 4 



(A9) 



Substituting Sx a obtained from Eq. (A7) and At from Eq. (A9) into a{x,x'), Eq. (3.101, gives a(x,x') as a function 
of Ax a , x a and the t„'s: 



<?(x,x') =- 



(2AJ - f) At 2 f (Ar - 2r Ti) Ar 



r — 2m 



(A6 2 + Acj) 2 ) - (2LA<j> + n- 2EAt) n 



1 f 2 / EM At 
2]j \f-2M ~ 

M (f Ti + Ar) A< 2 



LA<p) tiA 



M (r n - Ar) Ar 2 
r + -k +f 



(f - 2M y 

2(LA0- EAt + Tl )T 2 \ 



(ran + Ar) (A6 2 + A0 2 ) - 



2f T 2 Ar 
f-2M 



(A10) 
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If we now specify that x a coincides with the point where the world line intersects with the light cone of x, we can 
use the equation <r(x, x') = to solve for the t„'s in terms of Ai° and x a . This gives us 

rrn Ar 

ri =EAt - LA<fi + — ± p, 
2M — r 

±1 / {L 2 [Ar 2 + AM 2 (A9 2 + 3A0 2 )] - ALM 2 (EAt ± 2p)A<t> - 2EM 2 (3EAt ± 2p)At} Ar 



r-2 



i \ M 2 r 

ALM 2 f At 2 A(j) - AM 2 [f {EAt ± p) + 2Ar) At 2 + L 2 Ar 3 

Mf 2 

[ALM 2 (EAt - 2f Ar) A0 + 8E 2 M 3 f {AO 2 + Acj) 2 ) + L 2 Ar 2 - 2EM 2 (3EAt ± 2p)At] Ar 

M 2 (2M - f) 

{ (f ~ 2M) [ALM 2 r A(t) + AM 2 (Ef At =F r p + E 2 Ar) - L 2 Ar] + 8E 2 M 3 Ar} Ar 2 AL 2 MArAt 2 



M{f-2M) 6 

AT [AO 2 + Acj) 2 ) [{E 2 - 2) Ar - r (EAt - LAcj) ± p)] 



(All) 



with the higher order terms following in the same manner. 



A2(x,x') and cr a > (Eqs. (A7), (A8), (|3.16[) and (3.14) respectively), in terms of Ax a and x a . 



Using our equations for Ar, Eqs. (|A9[) and (All), we rewrite x a (and consequently 6x a ), the four velocity, u a , 



Bivector of Parallel Transport 

To calculate the bivector of parallel transport, g a b>(x, x'), we first write it in terms of a coordinate expansion about 



x, 



g a b '{x,x / ) = 5 a b , + G a bc {x)6x c ' + G a bcd (x)8x c ' 8x d ' + G a bcde (x)Sx c ' 6x d ' Sx e ' +..., 



(A12) 



where the coefficients G a b ...(x) are functions of x a written in terms of Ax a and x a . Calculating g a b '. c > [x, x') is straight 
forward: 



5 Vc {x, x 1 ) = G a bc (x) + 2G a bcd {x)Sx d ' + 3G a bcde (x)Sx d '5x e ' + AG a bcde {x)Sx d ' 5x e ' ' Sx? + ■■■ . 



(A13) 



Using the identity g a b >. c 'cr c ' = g a b >, c >cr c ' - ^w cl g d ^ c ' = with Eqs. ( |3.12[ ), ( |A13[ ), ( |A12| , and our expression for Sx a ' 
(obtained from the previous section), one can calculate the above coefficients and hence obtain the bivector of parallel 
transport, g a b '(x,x'), in terms of Ax a and x a . 



Scalar Singular Field 



Combining Eqs. (2.13) and (3.1), the scalar singular field can be written as 



$(s)( a 



A2(x,x') 

a c i (x, x')u c ' (x') 



X =I(,dv) 



T(adv) 



V{x,x(r'))dT. 



(A14) 



We already have everything required for the first term here, which gives the direct part of the scalar singular field. 
It should be noted that x' = X( ret ) and x' = X( a dv) are the equivalent of setting ±p — —p and ±p = +p respectively 
when substituting the r„, Eq. (All), into A 1 / 2 , a a i and u a . 



In the scalar case, Eq. (3.5) for the scalar tail part becomes 

oo 

V(x, x') = V n (x, x')<r n (x, x'). 



(A15) 
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To calculate coordinate expansions of the V n , first we require a coordinate expansion for Vq about x of the form, 

Vq(x, x) — vq(x) + voa(x)Sx a + VQ a b (x)Sx a 5x b + vo a b c (x)Sx a Sx b Sx c +••• . (A16) 



The 'initial condition' described by Eq. (3.6b) in the scalar case then becomes 

1 



2cr ;a V 0:Q ' - 2V A"2 C r a 'A2. a , + 2V a + (□' - m 2 - £R) A 2 = 0, 



(A17) 



and from this, it is quite simple to read off expressions for the coefficients vo a ... ■ Once we have Vq to the desired order, 
we compute a coordinate expansion for V n (n > 0) of the form 



V n (x, x') = v n (x) + v na (x)Sx a + v nab {x)8x a 5x b + v nabc (x)Sx a Sx b 5x c + 



The recursion relation for V„ , Eq. (3.6a), in the scalar case is then 

l 



2na' a 'v n -y - 2nV n A" 2 CT ;a ' A 2 ;a , + 2n {n + 1) V n + (□' - m 2 - £R) V n - X = 0, 



(A18) 



(A19) 



from which we can obtain expressions for the coefficients, v na . Here, the number of terms which must be computed 
is determined by the accuracy to which we require the singular field. For the present calculation, we require up to 

VQabcde, Viabc and V 2 a ■ 

Once we have V n to the required n, using Eqs. (3.10) and (A15), along with our expression for 5x a obtained from 
Eq. (A7), we get V(x,x') in terms of Ar, Ax a and x a . This can be easily integrated over r as required by Eq. (A14). 
Our final expression for $( s )(x) is then obtained by using Eqs. (A9) and (All) to remove the At dependence. As 
before T( rc t) and T( a d v ) are obtained by allowing ±p = —p and ±p = +p, respectively. 



Electromagnetic Singular Field 



For the electromagnetic singular field, we use Eqs. (2.18) and (3.1) to give 



A2(x,x / )g aal (x,x / )u a ' (x') 

<7 C < (X, X')u c ' (x') 



-| x =Z (adv ) 



T(adv) 



V aa ,(x, Z (T))u a '(x')dT, 



where V aa > (x, z(t')) is given by Eq. (3.5) 



V aa \x,x>) = Y,V™' (x,x')a n {x,x') 



(A20) 



(A21) 



n=0 



and the relevant metrics at x and x' can be used to lower indices. We require a coordinate expansion of V$ a of the 
form, 

V aa '(x,x') = vf{x) +vf b {x)8x b ' + vf bc {x)5x b ' 5x c ' + vf bcd {x)5x b ' 5x c ' 8x d ' +■■■ . (A22) 
Substituting this into the "initial condition" in Eq. (3.6b), which in the electromagnetic case is 

2a' b ' V Q aa ' , b> - 2V aa 'A-2V b 'A2> + 2V Q aa ' + (s a ' b ,a' - R a ' ' b >) ( A^g ab ' ) = 0, (A23) 



the coefficients of Eq. (A22), v?j a b ..., can easily be recursively obtained. It should be noted that the covariant 



derivatives do require the appropriate Christoffel symbols, which can be obtained from the suitable metric at x' . 
Next, we construct coordinate expansions for the V£ a ■ These have the form 



V™'{x,x') = <°» + b (x)5x b ' + bc (x)5x b ' 5x c ' + vf bcd (x)Sx b ' Sx c ' Sx d ' + 



b' x„c' 



„b' X„c' X„d' 



Substituting Eq. (A24) into the recursion relation (3.6a), which for the electromagnetic case becomes 

l 



2na' b 'V° a > -2nlC a 'A"2 CT ^'A2 ;b , +2n{n + l)V^ + f <5 a b >O f — R a b < \ V^—i = 



rab' 



(A24) 



(A25) 



we can recursively solve for the coefficients of Eq. (A24), v® a Once we have V£ a to the required n, we carry out 
the same remaining steps as in the scalar case and use Eq. (A20) to calculate the electromagnetic singular field. 



3G 



Gravitational Singular Field 



In the gravitational case, Eqs. (2.23) and (3.1) give 



;(S) 



2/x 



A2(x,x')g a , {a g b)b ,(x,x')u a (x')u» {x') 
a c > (x, x')u c ' (x 1 ) 



+ 2 M / tidV> Vata'v (x, z(T'))u a ' {x')u b ' {x')dT, (A26) 



where V a ba'b'{%-> Z { T ')) is given by Eq. (3.51. For the gravitational case, this is 

oo 

V aba ' b '(x,x') = Y,V: ba ' b '(x,x')<T n (x,x'), 



(A27) 



n=0 



where the appropriate metric at x or x' can be use to lower indices. The coordinate expansion for Vq 1 " 1 b (x, x') is of 
the form 

V aba ' b '(x, x') = v% ba ' b '{x) + v% ba ' b ' c (x)Sx c ' + v* ba ' b ' cd (x)Sx c 'Sx d ' + v% ba ' b ' cde (x)5x c ' 5x d ' 5x e ' + ■■■ . (A28) 



We replace V r afca b (x, x') with Eq ( A28) in the "initial condition" described by Eq. (3.6b), which for the gravitational 
case is 

2 t r c V abQ ' 6 ' ' , d - 2V afca ' 6 'A-5 (r c 'Ai c , + 2V aba ' b ' + (d a ' d 8 b> d ,W + 2C a ' c , b ' d ^j U^g c '^ a g^ d '\ = 0. (A29) 



This equation may be used to recursively solve for the coefficients of Eq. (A28), VQ ba b c .... Next, the coordinate 
expansion of V£ ba b (x, x') for n > has the form, 



VZ ba ' b ' [ Xl x') =v a n ba ' b ' (x) + v a n ba ' b ' c (x)6x c ' + v: ba ' b ' cd (x)Sx c ' 6x d ' + vf a ' b ' cde (x)Sx c ' Sx d ' Sx e ' + 



aba'b' 



aba'b' 



aba'b' 



c' i n aba'b' 



„c' X,„d' re' 



(A30) 



Substituting this into the recursion relation of Eq. (3.6a), which for the gravitational case has the form 



2na' c K 



• ' T -«'*//.' -),.T -aba'b' x ■.->„:.-' \ ., ■ >,,,, i > \r»im"b' ( \a' xb' _n' i nna' ^b' \ yabc^d' _ g (A31) 



2nV£ ba b A-2a' c A 2 ;c , + 2n (n + 1) V™ b + [S a c ,S b d ,W + 2C 



we can recursively solve for the coefficiens of Eq. (A30), v^ ba b c .... As in the previous two cases, once we have 
V£ ba b (x,x') for the required n, it is straightforward to calculate the singular field using Eq. (A26). 



Appendix B: Covariant bitensor expansions 



In this Appendix, we give covariant expansions for the beitensors appearing in the formal expression for the singular 
field, Eq. (2.6 ). These are given in terms of the biscalars s = (g ab + u a u b )aaCri, (the projection of a a (x, x) orthogonal 



to the worldline), and r = <7 a u° (the projection of <7a(x,x) along the worldline). In writing the coefficients, we use 
the notation [T ai ... an ](k) to denote the term of order e k in the expansion of the tensor T ai ... an , so that 



T 

ai •••a rl 



oo 

E 

fc=0 



[Ta 1 -a n \(k)( k ■ 



(Bl) 



Advanced and retarded points 



Eq. (2.6) for the singular field includes bitensors at points x' on the world-line between the advanced and retarded 



points of x. We consolidate this dependance to a single arbitrary point, x, on the world by expanding the dependence 
on x' about x. Denoting the proper distance along the world-line between a^advj/^rct) aud x by Ar, we may write 
the expansion of this distance in powers of e as 

A _ | _ . „ A if i s)" 2 R U(Jua if i s) ((r i s)R U aua;u — R U crua;a) 

Ar (1) =r±s, Ar (2 ) = 0, Ar (3) - T ^ , At w = T ^_ , 
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Ar (5) = T 



(f ±s) 2 



360s 3 



^ ^RuauaRucrua (? ~F 3s)(r zb s) -f- S (r zb s) ^3-R ucrucr;liU -h ^RuauaRucru 



At, 



(6) 



± 



(r ± s) : 
4320s 3 

i t2 



-(f± s)(s 
|30i? MC 



}■ 



^Ruaua\ua ^^RucruaRua a J H~ ^Ruaua;acr ~b ^Rua a aRua a 
Ruaua-air T 3s)(f ± s) - iJ MM; „(f =p 4s) (f ± s) 2 
^Ruaua^aaa ^^Ruaad-.a Ruaa 
(v zb 5) ^3/? W(7Ui(7 ; uerer 3QR u &{ja :( j R U(7U ~ Ruaa (]-6Ru(7ua:(j H~ ^Ruaua-.a 10-^iitrcrd;i 



H~ 2(r zb 5) {^Ruaua\uua ^^Ru(Tua,uRuacr ~b Ruau i^^Ruaua-.a ~\~ 
6(r zb s) (yRu(7U(7\uuu ~b ^Ruaua\uRuau ) }■ 



(B2) 



Advanced and retarded distances 

Taking two derivatives of the world function, we obtain a bitensor that has the covariant expansion 

^ab = <?afe — \Racbd& C <J d + j^Racbd.e^ V° ^ + (^RacpdR? ebf + ^Racbd:ef^j^ C ^ d G e G^ 
^{^Y^Racpd^eR^ fbg ~b ~y2£)RacpdR? 'ebf;g H~ ~^^Racbd,ef g^j^ ® (J (J® 
(^945 ^acpdR^ eqf R^ gbh H~ T^Q^Racpd-.ef R^ gbh ~b 5040 ^ ac pd>;eR^ fbg:h 
+ r^lRacpdR P ebf-gh + ^^Racbd-ef gh^j^ &° ^ ' G^ G^ ' G h 

~b ^ 20160 ^ ac pd\eR^ fqg^^ ' hbi 302^ RacpdR^ eqf ;gR^ hbi H~ 30240 -^&cpd;efgR^ hbi 

~b 20160 ^ ac pd;ef R^ gbh;i ~b 20160 ^ ac pdR P eqf R^ gbh;i ~b 20160 ^ ac pd;eR P fbg;hi 

I 11 p pp 1 1 p \ —C d e f g h i 

30240 ac pd e bf;ghi ~r 20160 o,cbd;efghi J G G G G G G G . 



(B3) 



For the singular field, we require the expansion of [cr a <it a 1(z±, x ). Writing [<T a m a ](r) = [<7 a /it a ](z(t),x), expanding 
the dependence on r about 5 (using the method of Sec. Ill A and making use of the above expansion of a ab ) and 



evaluating at r = r±, we obtain the coefficients of the expansion of [a a 'U a ]± = [a a 'U a ]{z±, x) about a;. They are: 



r (i) 



r (3) = 



Gs 



FL, 



"(4) 



r ± s 
24s 



17 



■ucruacr 



(r =F s) - R u(7uau (r ± s)(r =F 2s) 
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"(6) 



360s 3 



q_2 ^ ) ^Ruauaacr ~t~ ^Ruaa Ruacra) ij" i (^* ~F 2s) (S-R^fj^fjer^ 16-/? uo -(7 Ruaua) 

-I 2 

5 (f 2 - s 2 )i? M(TMCT j, 

6(r zb 5) (V ^- 4s) (^Ruauauuu H~ ^RuauduRuau 



L UO"U<TUU 



Hl<ju<j<t<7(j 



AR a R -1 

QR U n-rr (jRl, 



+ (f±s) 3 (f=F3s)(3i?. 

2(r zb s) (r ^ 2s) ^S/^ct^ct^cj-ct ^^RuoadcrRuau Ruaa {^^Rucruaa "T~ ^Ruauad ~~ 10-R 

^ORu&uauRucra H~ Ruau iX^Ruauda ~b ^Ruauad ~~ ^^Ruaadu) 



■2(f±s) 3 (fT3s) 



+30i? 



(r 2 - s 2 ) 2 R uauaa - (f ± s) 3 (r 2 =p 3rs + 4s 2 ) J R U(TmT „ |. 



(B4) 



Van Vleck Determinant 



The Van Vleck determinant has the covariant expansion 

A 1 / 2 ^^') = 1 + ^R ab a a a b - ^R ab -ca a a b a c + ^R paqb RP c i d + ^RabRcd + ^Rab-,cd)«" 
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{^r^RpaqbR P c Q eke + ^^R a bR C d;e + 7^-Rab;ede) a a (J b (J C (J d a C 
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Writing A 1 / 2 (r) = A 1 / 2 (z(r), x), expanding the dependence on r about x (using the method of Sec. Ill A) and 
evaluating at r = t±, we obtain the coefficients in the expansion of A^/ 2 = A 1 / 2 (z±,a;) about x. Specialized to the 
vacuum case, they are: 



A 



1/2 _ 
(o) _ 



1, A}/ 2 = 0, A^ = 0, A 



1/2 _ 
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Bivector of parallel transport 

The derivative of the bivector of parallel transport has the covariant expansion 

g a a g a 'b:c{x, X 1 ) — —-RbacdP d + QRbacd;e^ d (J e — — (^RbapdR P ecf + Rbacd;efj O d a e U i 
+ {^ ) RbapdR P ecf -g + 7^Rbapd;eR P f eg + ^Rbacd-.e f g^ <J d G S O 9 



(^2^0^ bapd ^ P ec f' 9h 120^ bapd,e ^ P ^ cg,h igo^ bapd ' e f ^ ' gch 
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"240 



RbapdR P eqfR q gch + Rbacd;efgh ) CF d C (J f (T 3 (J h 



(B7) 



For the singular field, we require the expansion of g aa ' U a (z± ,x). Writing [g aa 'U a ]( r ) = [gaa' u<1 ](z(t),x), expanding 
the dependence on r about x (using the method of Sec. Ill A and making use of the above expansion of the bivector of 



parallel transport) and evaluating at r = t±, we obtain the coefficients of the expansion of [g aa 'U a ]± = g a a'U a (z±,x) 
about x. They are: 
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(4) 24s 



^ucrua-.uu 



Ruoua\aa ~b RuaabRucra "f (f i s) Ru 
+ (f zb s) \ ' i^R^aab ^Rucrab) ~ Rucrua-ua + Rudu {R U crab + (f" i s)R uau l)^j J J, 



zb- 



1 



(5) 2160s 



ffa s (f±s-){s[l8( 



R 



ui7ua\aai7 ±L uaua\ua{j 



Ri, 



r {f±s) + R 



uaua-.uua 



(f±s) 2 -R 



ucruatuuu 
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~\~QRuaab {^J Ruaa ;tr H~ if i s)(3i?- uau -(j ^Ruaa :u H~ ^Ruaua^ ) QRuau \u(j' i s) ^ 

H~9-^uoorj[> ^4J2 ucr(T ^ H- (r zb s s ){hR uau 3R U(7( j - u ) — QRucru \uir *"0 ^ 
+2i? (21(f ± s)(2i?„ aCT b ;ff + 3i? UCTa & ;CT ) 

+ l6R uaua > ) - 21i? UQU ;u (r ± s) 

+54(f ± s)i? U( xab (^Ruau -.a ~ 1R uau h , u (f ± 
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3R 
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Scalar tail 



The scalar tail bitensor, V^(:e,:e'), may be expanded in a covariant series by writing it in the form of a Hadamard 
series, 

V(x,a/) = F (x,a; / )+Fi(a;,a; / )o-(a;,a; / ) + -- - , (B9) 
and expanding each of the Hadamard coefficients Vo(x, x'), Vi(a;, a;'), • • • in a covariant Taylor series, 

Vo = «o ~ \vq. c u c + \vQ cd a c u d + \ (-|wo(cd ;e ) + iWo ; (cde)) cr c o-'V H , 

Vi = ui - \vy- c o c -\ . 

The series coefficients required to obtain the expansion of the singular field to C(e 4 ) [V(x, x') to C(e 3 )] are given by 
y o = |((€- |)i? + m 2 ) , 

For the singular field, we require the expansion of f T(adv) Vdr'. Writing V^t) = y(z(r),x) and expanding the 

T 

III A| and making use of the above expansion of V), we obtain an 
expansion in powers of At that can be trivially integrated between r = t_ and r = t + . Specialized to the vacuum 
case, the required expansion coefficients are then: 



T-(adv 



(i) 



0. 



T(adv) 



Vdr' 



(2) 



0, 



T(adv) 



'(ret) 



(3) 



0, 



T(adv) 



(4) 



0. 



Electromagnetic tail 

The electromagnetic tail bitensor, V a b'(x, x'), may be expanded in a covariant series by writing it in the form of a 
Hadamard series, 

V aV [x,x') = gv b [V 0ah (x,x') + V lab (x,x')a(x,x') + •••], (BIO) 
and expanding each of the Hadamard coefficients, Vo a i,(x, x') 1 V\ a b{x, x'), • • • , in a covariant Taylor series, 

V ab = Vo( a b) + (-^W0(a6);c + «0 [ ab ] c ) <J C + \ (f (ab)cd ~ Va[ab](c-d)) ^ '<* 

+ \ (~| W 0(a6)(cd;e) + \vq (ab):{cde) +V Q [ a b]cde) ^ H , 

Vi a b = Vi ( afc ) + (-|wi(a6);c + V\ [ab]c) 0" C + • • • . 
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The series coefficients required to obtain the expansion of the singular field to 0(e 4 ) [V a b' (x, x') to C(e 3 )] are given 
by 



v 0(ab) — \Rab ~ J^Rgab, 
Vq [abf 
Vq (ab) 



Vq [ab] 



cde 



6 RC [b;a], 

1 p ;(ed) | 1_ p ped j 1_ p pqc p d 

Q^ab i 12 a b ' 12 ( a ^fypq 

I I _ 1 p c rtpqrd 1 pc d ppg , 1_ pc pdp J_ p pcd J_ p:cd 

i yab \ iso PQ r 180 P 9 90 P 7 2 - n '- rt 40 

d;e) J_ p p(c p <je)g , J_ p (c;d pe)p 
3Q Jt ao ;g Jt p ' go Q bp ±tj 



cd\ 



-i? (c r M de) - — i? (c r uR de) - —Rs pq( - c Ru 
20 1 [a;b\ yi [a;bj-"- 20 l a -"-bjpg 



t ( c R de )<p Lb ,p( c /? d ; 

>0 Jt aop JL 20 a b ±h p 
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and 



^1 (ab) — ^RapqrRb^ 4~ gRapRb^ 2A RRj(1 b 24 ^-^ab 



■9ab ( 720 RpqrsR 1 " 1 720^ o ^ P<? + 288^ ^~ 12o'~'^) ' 



^1 [ab] 



1 p par p 
t -i n -*Ma "b 



240~ L 1 a -^bjpgr i 24 [a -^bjpjg ' 120 [ a b]q;p 

T2Q R pq[aR^^ M 2A R ^ ;[a-^b]p ~t~ 24-^^a-^b] ;p 360^^ -^p[a;b] 24 R^ [o.Rb]p' ~H j^RR [a;b] 

+ _J_ R c r P ?-/?,, RWC ?-/?,, B r <3P c 1_ D c ppg , _J_ p pep;g 1_ p pen 
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1 on ^Ma -'•■blarr. 



These are the same as those given by Brown and Ottewill [Blj with the exception of V\ [ a b] c , where we have corrected 

a sign error in one of their terms and combined another two terms into a single term. 

For the singular field, we require the expansion of J" r < adv) V a b'U b dr' . Writing [V a b'U b ]{t) = Wab'U b ](z(t),x) and 

( rct ) 

expanding the dependence on r about x (using the method of Sec. Ill A and making use of the above expansion of 
V a b' ) , we obtain an expansion in powers of At that can be trivially integrated between r = r_ and r = t + . Specialized 
to the vacuum case, the required expansion coefficients are then: 
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Gravitational tail 



The gravitational tail bitensor, V aa >bb'(x, x '), mav be expanded in a covariant series by writing it in the form of a 
Hadamard series, 

V aa 'bb'{x,x') = g a ' C gb' d [V acbd(x,x') + Vi acbd (x, x')u(x, x') H ], (B12) 

and expanding each of the Hadamard coefficients, Voacbd(x, x'), Vi ac bd(x, x'), • • • , in a covariant Taylor series, 



VqAB — v (AB) + (— 2 V 0(AB);e + v 0[AB]e) ^ + \ ( v 0(AB)ef ~ v 0[AB}(e;f)) & 

+ I (-^ V 0(AB)(ef;g) + J«Q (AB)-(efg) + VQ[AB]efg) <J e (J S O 9 -\ , 

V\AB = Vi (AB) + \—2 V l{AB);c + V\ [AB]e) ^ H • 



The required coefficients for Vqab are 

VQ(AB) = V Q(aJ^d) = "Cacbd, 
Vo[ABf = 0, 



(B13) 
(B14) 



Vq (AB) 



ef _ „ ef _ ,.;(ef) _ \r< P(e/-i f) i 1_ nj>l{ e n, f) _ _J_TT , ,a e ffVW*C fRI ^ 

— V 0(abcd) ~ q^acbd ^<^ac ^bdp T ^Qac^b ^ dpq Tf^ ll a.bcd9 ^ ^pqrs, (OLO) 



720 



and 



where 



«0 [AB] 



e/9 — 



W 0[a6cd]" J3 — Y^5acO[ b '" 1 ^ ,J ^ d^q"' + -^9ac^bdpe;qC P f q g , 



1 1 

Ilabcd — ^9ac9bd + ^9ad9bc + K 9ab9cd- 



(B16) 



(B17) 



In Eqs. (B13) - (B16), the right hand sides are understood to be symmetrized on the index pairs (a&) and (cd). The 
required coefficients for Vi ^ are 

«1 (AB) = «i (^fc^d) = T^nCacftd + -C a P b q C cp dq + ^C ac Pq Cbdpq - 

and 



0(j 



gacgbd^'C. 



pqrs 



l 
720 



pqrs 



(B18) 



— ^2 a P b q C C pdq' e — C P d q C a pbq' e ) + — (C a Peq Cbcdp;q ~ C ' c peq 'C 'dabp;q) + ^SacC"^ 9 ' Cbdpq\ri (^19) 



where again there is implicit symmetrization on the index pairs (ab) and (cd). When k = —1/2, Eq. (B18) agrees 
with Eq. (A23) of Allen, Folacci and Ottewill [88] specialized to the vacuum case. Our expressions also agree with 
Anderson, Flanagan and Ottewill 89J, but we write them here in a slightly more compact form. Note that the 



expressions (B13) - (B16) and (B18) - (B19) are all traceless on the index pair (cd), aside from the terms involving 
the tensor H a bcd- This means that performing a trace reversal on the index pair (cd) is equivalent to changing the 
value of k from to —1/2. For the calculation of the gravitational singular field, we require an expansion of the trace 
reversed singular field and so we choose k = 0. 
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